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This article gives a comprehensive review on thick brane solutions and related topics. Such models 
have attracted much attention from many aspects since the birth of the brane world scenario. In 
many works, it has been usually assumed that a brane is an infinitely thin object; however in more 
general situations, one can no longer assume this. It is also widely considered that more fundamental 
theories such as string theory would have a minimal length scale. Many multidimensional field 
theories coupled to gravitation have exact solutions of gravitating topological defects, which can 
represent our brane world. The inclusion of brane thickness can realize a variety of possible brane 
world models. Given our understanding, the known solutions can be classified into topologically 
non-trivial solutions and trivial ones. The former class contains solutions of a single scalar (domain 
walls), multi-scalar, gauge-Higgs (vortices), Weyl gravity and so on. As an example of the latter class, 
we consider solutions of two interacting scalar fields. Approaches to obtain cosmological equations 
in the thick brane world are reviewed. Solutions with spatially extended branes (S-branes) and 
those with an extra time-like direction are also discussed. 
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I. INTRODUCTION 

At present, the presence of extra dimensions is an integral part of almost all physical theories playing a fundamental 
role in trying to explain physical interactions based on common principles. Also, any realistic candidate for a grand 
unified theory, such as superstring/M-theory, should be multidimensional by necessity, otherwise, it will contain 
undesirable physical consequences. Traditionally, it has been considered that in such theories the observable 4- 
dimensional spacetime appears as a result of compactification of extra dimensions, where the characteristic size of 
extra dimensions becomes much less than that of 4-dimensional spacetime. But, recently, a new idea, called brane 
world, has come up, where we live on a thin leaf (brane) embedded into some multidimensional space (bulk). Because 
of the lack of experimental and observational data, there is no reason to prefer a particular class of multidimensional 
models of gravity. Actually, various kinds of multidimensional theory have been studied in modern theoretical physics. 

In this article, we will review brane world models where our brane universe has a finite thickness^. The reason is 
mainly two-fold: Firstly, such solutions appear in various multidimensional field theories coupled to gravity, leading 
to a variety of possibilities of brane world, which are interesting by itself. Secondly, brane thickness should be an 
essential ingredient to extend the idea of brane world to multidimensional spacetimes. 

Before coming to our main subject, it is important to keep the main history of multidimensional gravity in mind. 
Thus, we briefly discuss the history of multidimensional gravity and especially why the idea of brane world arose. 

A. Prebrane epoch of multidimensional theory 

1. Why is the visible universe 4- dimensional? 

After the birth of Einstein's general relativity, the natural question about the dimensionality of the world which we 
live in appeared. Within the framework of Einstein's gravitational theory, space and time are unified and it allows us 



^ In this review, we will basically follow the notations, sign conventions and definitions of fields which were adopted in the original papers. 
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to realize that the surrounding world is a 4-dimcnsional one. This statement is a consequence of a combination of the 
principle of Occam's razor and Einstein's idea about the geometrical nature of the gravitational field. Hereupon, the 
gravitational field is not an external field on the background of a specified spacetime, but a geometrical characteristic, 
i.e., metric, self-controlled on any sufficiently smooth manifold. 

One of the questions inevitably arising in such a case is about the dimensionality of our world. Naturally, it is 
necessary to consider two questions concerning the dimensionality of the physical world: 

• Why the physically observed dimensions of our Universe = 3+1 (space + time)? 

• If the real dimensions are still more than 4, why extra dimensions are not observable, and what about any 
consequences of multidimcnsionality for a 4-dimcnsional observer? 

Strictly speaking, the answer to the second question could serve as an answer to the first question. Indeed, if extra 
dimensions are hidden from direct observation then our world looks 4-dimensional effectively. 

Although at the present time, it is not possible to give an exact answer to the first question about reasons of the 
4-dimensionality of our world, investigations in this direction show that the 4-dimensional spacetime is marked out 
somehow among spacetimes with a different number of dimensions. Let us enumerate these properties: 

• Einstein's equations in Lorentzian manifolds with dimensions d < 4 lead to a flat metric. 

• Maxwell's equations are conformally invariant only in a four-dimensional spacetime. 

• In Newtonian gravity, circular orbits of point-like bodies in cental-symmetric fields are stable only on manifolds 
with dimensions d < 4. 

• In spacetime with dimensions d > 4 the Schrodinger equation either does not have any energy level, or an energy 
spectrum is unbounded below. 

• Huygens principle is valid only for spaces with odd dimensions. 

• Quantum electrodynamics is rcnormalizable only in a spacetime with dimensions d < 4. 



2. Kaluza-Klein theory 

The present stage of the development of multidimensional theories of gravity began with Kaluza's paper [l[ where the 
foundations of modern multidimensional gravitational theories have been laid. The essence of Kaluza's idea consists 
in the proposal that the 5-dimensional Kaluza-Klein gravitation is equivalent to 4-dimensional Einstein gravitation 
coupled to Maxwell's electromagnetism. In the primordial paper there was the proposal that the 655 metric component 
is a constant, and it is not varied when obtaining the gravitational equations. Later on, this proposal was excluded. Let 
us show the basic advantages and disadvantages of the standard 5-dimensional Kaluza-Klein theory. The advantages 
are: 

• The 5-dimensional Einstein equations could be reduced to a form identical to the equations of 4-dimensional 
gravitation interacting with electromagnetic and scalar fields. 

• The geodesies equations in a 5-dimensional space lead to the equations of motion for charged particles in the 
electromagnetic and gravitational fields. 

• If one supposes that all quantities in the Kaluza-Klein theory do not depend on the 5-th coordinate then it 
transforms as a transformation of coordinates in the following way: 

a;'' = .T^ + /(.t"), M = 0,1,2,3 (1) 
This coordinate transformation leads to the usual gradient transformations in Maxwell's electromagnetism. 
The disadvantages are the following: 

• Why is the 5-th coordinate unobservable in our world? 

• Einstein noted [Ij that it is not possible to give a natural physical meaning to the metric component G55. 

• The equation for a scalar field 7^.55 — G^^TZ = gives the stiff and unusual connection between the 4-dimensional 
scalar curvature R and the invariant Ffj^^F^^" of the electromagnetic field. 
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3. Multidimensional theories 



Next, multidimensional Kaluza-Klein theories were considered. One of the most important questions is the question 
about the unobservability and independence of 4-dimensional quantities on extra dimensions. The most general 
approach for a solution of this problem consists in obtaining the solutions of the multidimensional Einstein equations, 
which yield characteristic linear 4-dimensional spacetime sizes of many orders of magnitude bigger than linear sizes 
in extra coordinates. It is also necessary that the 4-dimensional metric of the obtained equations does not depend on 
extra coordinates. This could be done by refusing vacuum equations, and introducing torsion or quadratic terms 
of curvature Q, or by adding multidimensional matter. The last variant called "spontaneous compactification" of the 
extra dimensions was suggested by Cremmer & ScherkJ^Q. Unfortunately, this approach gave up Einstein's original 
idea about the geometrization of phy sics. Following [10|, multidimensional Kaluza-Klein theories can be classified 
into compactified theories [sj-fol IllL [l2j and projective theories [T3l|-(l5l|. 

According to the definition, compactified theories should satisfy the following properties: 

• some periodicity condition with respect to extra coordinates is fulfilled; 

• extra dimensions form some compact manifold M; 

• the linear sizes of this manifold M are small in comparison with the sizes of spacetime [ll|, [l^l . 

It is necessary to note that the metric which satisfies the above properties is a solution of the vacuum (4-|-d)-dimensional 
Einstein equations only if M is a flat manifold [3 • Manifold M with constant curvature can be a solution of the 
Einstein equations only in the existence of multidimensional matter torsion or quadratic terms of curvature 

0. In the projective geometry, a point is a set of projective points for which all relations of (n -I- l)-homogeneous 
coordinates are the same. In such a case, each projective line is a point in a real space. Such an approach leads to 
the conclusion that the metric does not depend on an extra coordinate. Projective field theories do not agree with 
experiment. At the present time, they are not considered as realistic candidates for an unified field theory. 



4- Problem of unobservabiUty of extra dimensions 

In the last few decades, the idea that our space, indeed, has a number of dimensions more than four became quite 
popular again. It is supposed that many problems of elementary particles physics could be solved by introducing 
strings, supersymmetry, and by increasing the number of spacetime dimensions. In increasing the number of di- 
mensions, a legitimate question arises: why physical properties of extra dimensions differ strongly from those of the 
observed 4-dimensional spacetime. Two solutions to this question were suggested, namely a mechanism of spontaneous 
compactification and brane world models. 

In the case of compactified extra dimensions, it is supposed that these dimensions are rather small and unobservable. 
In the papers (T7j-[23j the modern approach to the problem of the unobservability of compactified extra dimensions 
was introduced. In a d— dimensional spacetime, it is possible to consider gravity coupled to matter. For the purpose 
of compactification, a special solution of the block form M'^ = x Af''"^ is searched for, where M'^ is the 4- 
dimcnsional spacetime and M'^~'^ is a compact space of extra dimensions. and M'^~'^ are spaces of constant 
curvature. The block representation of the metric is consistent with the gravitational equations only if one supposes 
that the energy momentum tensor takes the block form T^ii^ ~ ~7i5/^iy in M^, and Tm„ = 72gmn in M'^~'^, where g^^^ 
and gmn are metrics of and M'^~^. There are various mechanisms of spontaneous compactification: Freund- Rubin 



compactification [2J| with a special ansatz for antisymmetric tensors; Englert compactification [25| setting a gauge 
field in an internal space equal to the spin connection; suitable embedding in a gauge group [26| ; monopole or instanton 
mechanism [23]; compactification using scalar chiral fields [1^, and compactification using radiative corrections [29| . 



B. Brane world scenario 



In this subsection, the new approach to the problem of unobservability of extra dimensions, called brane world 
scenario, is explained. This approach is quite different from the traditional compactification approach and allows 
even noncompact extra dimensions. According to the new approach, particles corresponding to electromagnetic, 
weak and strong interactions are confined on some hypersurface (called a brane) which, in turn, is embedded in 
some multidimensional space (called a bulk). Only gravitation and some exotic matter (e.g., the dilaton field) could 
propagate in the bulk. It is supposed that our Universe is such a brane-like object. This idea for the first time was 
formulated phenomenologically in f30j~[33|, and got confirmation within the framework of string theory. Within the 



5 



brane world scenario, restriction on sizes of extra dimensions becomes weaker. It happens because particles belonging 
to the Standard Model propagate only in three space dimensions. But Newton's law of gravitation is sensible to 
existence of extra dimensions. 

The prototype idea of brane world appeared a rather long time ago. Now the term "brane model" means different 
ways for solving some fundamental problems in high-energy physics. The pioneering works in this direction were done 
in the papers [30| and [3]| . These prototype models have been constructed phcnomcnologically and considered as 
topological defects, i.e., thick branes in the modern terminology. 

In [30|, an idea that extra dimensions could be noncompactified was suggested. It was supposed that we live in 
a vortex enveloping multidimensional space. Akama showed that "... The Einstein equation is induced just as in 
Sakharov's pregeometry . . . " . The main idea of this paper consists in using the Higgs lagrangian in a 6-dimensional 
flat spacetime 

C = + DM<P^D^'(j> + a|0|2 - b\cf>\' + c, (2) 

where Fmn = Qm-^n — QnAm and Dm4' ~ + ieAi\i. It allows one to obtain the induced Einstein gravitation on 
a vortex. The corresponding field equations have the vortex solution 

9 9 

Am = eoi23MArA(r) — ,0 = ip{r)e"'',r' = {x^Y + {x^ , (3) 

where A{r) and (p{r) are the solutions of the Nielsen-Olsen differential equations describing the vortex. Our Universe 
(the vortex) is localized inside the region 0{e){e = l/y/a). Further, introducing curvilinear coordinates 

Xm = Y^^ (x'') + n^x"\ i\/ = 0,1,2,3,5,6, 0,1,2,3, m = 5,6, (4) 

where X^ are the Cartesian coordinates and are the normal vectors of the vortex, one can show that the Einstein 
action is induced through vacuum polarizations. 

In the paper [3l| the idea that we live on a domain wall (brane in modern language) and "... ordinary particles 
are confined inside a potential wall . . . " , was suggested. The quantum toy model with Lagrangian 

£ = isA'pS^^ - imV - ^V'", A = 0,1,..., 4 (5) 

is under consideration, which describes a scalar field in the (4+l)-dimensional spacetime M'-^'^-' with metric qab = 
diag(l, — 1, — 1, — 1, — 1). The classical field equations for the scalar field ip admit a domain wall solution (which is a 
(l+l)-dimensional kink) 

This kink solution can be considered as the domain wall. In the paper the possibility of trapping particles with spins 
and 1/2 on the domain wall was also considered. 

For particles with spin 0, the linearized equation of motion for the field (p' ~ (p — pf^ is under consideration 

-aAa^¥''-"^V-3A(v?"0V' = 0. (7) 

It was shown that there exist three types of perturbations: 

1 . the first one is of particles confined inside the wall and with the energy E — k^; 

2. the second one is of particles confined inside the wall, but with the energy E = k'^ + |m^; 

3. the third one is not confined inside the wall and at large x'^ the energy is 

The possibility of interaction of the first type of particles with creation of the third type of particles and with 
subsequent escape of them beyond the domain wall was discussed. 

Furthermore, the question about trapping fermions on the domain wall was considered. The Lagrangian of massless 
fermions is 



(8) 
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It was shown that for such an interaction between a fcrmion ^' and a domain scalar field (y9, the Dirac equation has a 
solution 

where ip{x^,x) is a left-handed massless (3+l)-dimensional spinor. It was also shown that there exist two types of 
fermion perturbations: 

1. the first one is of fermions ([9]) confined inside the wall; 

2. the second one is not confined inside the wall and the energy exceeds hm/V^. It can be created only in high 
energy collisions. 

Note that in both papers [30| and [3l[ the idea about noncompactified extra dimensions was suggested. The difference 
consists in that in [sof our 4-dimensional world is the interior of a vortex, while in [sij our 4-dimensional world is a 
domain wall. 

We should also mention interesting earlier attempts in the context of supergravity and superstring theory. In [s^ , 
models of domain walls in A'^ = 1 supergravity theories were considered. Using this model, it was shown that there 
exists a class of domain wall solutions which need not to be Z2 symmetric. Such a solution describes a stable domain 
wall that divides two isolated but non-degenerate supersymmetric vacua, and at least one of them is an anti-de Sitter 
one. Moreover, it was shown that in supergravity theories gravity plays a non-trivial and crucial role for topological 
defects (in particular, for domain walls). 

Ref. [36| examined the possibility that some of the internal dimensions, at relatively low energy scales, may be of 
the order of a few TeV, within the framework of a perturbative string theory which could relate the supersymmetry 
breaking scale to the size of the internal dimensions. 

Although several ideas have been presented during the 80s and the early 90s, one of the most striking facts which 
activated the studies on brane models was the development in superstring theory and M-thcory since the mid of 90s, 
especially the discovery of D-brane solutions [13, HI] . A valuable contribution to development of this model was made 
in the papers [s^ and [i^ , following an earlier idea from [s^ . The authors considered a flat bulk geometry of (4 -I- d) 
dimensions, when d extra dimensions are compact with radius R. All the Standard Model particles are assumed 
to be localized on the branes at the boundary of the d-dimensional compact space. But the gravitational interaction 
feels the presence of extra dimensions on the distance scales less than R and differs from Newton's law there. By 
integrating over the internal space, it turns out that the 4-dimensional Planck mass Mp\ and the (4 -I- d)-dimensional 
Planck mass Mfund are connected by the following relation: 

^^pi - Ki.'X- (10) 

Since Newton's law has been checked up to the distances of the order of 0.1 mm, R can be of order 0.1 mm or less. 

In the papers [4l| - |43j the model of the Universe expanding as a 3-shell in a in 5-dimensional spacetime was considered. 
It was shown that the model can solve the hierarchy problem. Also the relation between the size of the Universe and 
its thickness was found. 

In the papers [i^l and [i^, Randall and Sundrum considered a non-flat bulk space. The metric has the form 

ds2 ^ e-2*"'<=l^l77^^dx'^da;'' - rjdy^ , (11) 

where fc is a parameter of order of the fundamental Planck mass; Tc is the radius of a 5-th coordinate. This metric is 
a solution of the 5-dimensional gravitational equations following from the Lagrangian 

S = ^f d^xdy^{-R-X)+ J2 j d-'x^-h^'^ - y^*)) , (12) 

i— 1,2 

where hfu, C^"^ and V^"^ are the 4-dimensional metric, the Lagrangian and the vacuum energy on (z)th-brane, re- 
spectively. One can show that the relation between the 4-dimensional Mpi and the fundamental A/* Planck masses 
is: 

Ml, = % [1 - e-^^-^-] ^ ^ . (13) 
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Namely, the scale of physical phenomena on the brane is fixed by the value of the warp factor. On the brane near 
y = TT (observed brane), the conformal factor for the 4-dimensional metric is = q-'^^^o-^ ^ q;^^ physical masses should 
be calibrated by taking this factor into account. For example, if kvc ^12 then masses of order TeV could be generated 
from the fundamental Planck mass ~ 10^^ TeV. 

Various physical problems have been considered in five-dimensional (and even multidimensional) brane models, 
especially in the Randall-Sundrum models and their extensions. See e.g.. [46l - l50j for reviews. In many works, it was 
assumed that the brane is infinitely thin. Although in the thin brane approximation many interesting results have 
been obtained, in some situations the effects of the brane thickness cannot be neglected. 

C. Thick branes 

From a realistic point of view, a brane should have a thickness. It is also widely considered that the most fundamental 
theory would have a minimal length scale. In some cases the effects of brane thickness can be important. The inclusion 
of brane thickness gives us new possibilities and new problems. In many multidimensional field theories coupled to 
gravity there are solutions of topological defects. They have lead to a richer variety of brane worlds (see sections HIl 
and IIII[) . At the same time, brane thickness also brings ambiguity in the definition of the effective four-dimensional 
quantities (see section W\ . 

We now give our precise definition of thick branes, to avoid possible problems related to possible differences in 
terminology. Our definition is based on the following form of a multidimensional metric: For five-dimensional problems, 
the solutions whose metric arc 

ds^ - a\y)gf,,dx^dx'' - dy^, (14) 

where — oo < y < oo \s the coordinate of the extra dimension, are considered. The four-dimensional gin, is Minkowski 
or de Sitter (or anti-de Sitter) spacetime. a{y) is the warp function, which is regular, has a peak at the brane and falls 
off rapidly away from the brane. A typical behavior of the warp function in thin and thick brane solutions is shown 
in Fig. [TJ When the Z2-symmctry is assumed, a{y) = a{—y). The normalizability of the graviton zero mode gives the 
condition that dya{y)'^ is nonvanishing and finite. In five-dimensional problems, the thin brane approximation is 



warp factor 




thin-brane 
thick-brane 



FIG. 1: The typical behavior of the warp factor a{y) in the 
thin and the thick brane solutions are shown. The dashed line 
denotes the warp factor of the thick brane, while the solid one 
does that of the thin brane. 

valid as long as the brane thickness cannot be resolved, in other words, if the energy scale of the brane thickness is 
much higher than those in the bulk and on the brane. In contrast, when thickness becomes as large as the scale of 
interest, its effect is no longer negligible. 

With more than six-dimensional spacetimes, ti > 1, it is assumed that the metric function has the form 

ds^ = a^{r)gi,^dxt'dx'' - b^{r){dr^ + r^dnl_i). (15) 

The four-dimensional metric g^^i, also represents Minkowski or de Sitter spacetime. dfl^_i is the line element of an 
(n— l)-dimensional compact manifold (the unit (n— l)-spliere in most cases). The extra coordinates are composed of 
the radial one r (0 < r < oo), and those of the {n — l)-dimensional ones. It is possible to put the brane at the center 
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r = by an appropriate coordinate choice. Here a(r) is the warp factor, which has a peak at the brane. The functions 
a(r) and b{r) are regular everywhere and the integral Jp dra'^b"r"~^ is finite, in order to obtain a localized graviton 
zero mode. In such higher-codimensional cases, the situation is quite different from that in the five-dimensional one. 
In the context of the thin brane approximation, in approaching the brane, usually self-interactions of bulk matter, 
e.g., gravitation, are divergent. Then, it is impossible to describe the matter interactions localized on the brane. To 
study such models, brane thickness becomes an essential ingredient which plays the role of an effective UV cut-off. 

Note that it is obvious that such a definition is rather ambiguous: for example, one can redefine the coordinate r 
in the following way 

^/h(r)dr = dz, (16) 

and then the brane metric can be rewritten as follows 

ds^ = a^{z)gf,^dx''dx'' - [dz^ + r^{z)dnl_^] . (17) 

In sections ITIV, exact solutions of static thick brane solutions in various classes of field theories will be considered. 
In sections V-VI, time-dependent situations are considered, and then the bulk metric is no longer the form of Eq. 
(El) or Eq. ([HI). 

D. Purpose and construction 

In this paper, we will give a comprehensive review of works devoted to thick brane solutions, and related topics. 
To our knowledge there has been no review on this subject although an enormous number of works have already been 
done. Thus, now is an appropriate time to collect these works. We hope that this paper becomes a good starting 
point for new studies. 

In our understanding, all known thick brane solutions can be classified into two large groups: static solutions 
(presented in Sections II- IV) and solutions depending on time (presented in Sections V-VI). In turn, the static solutions 
could be classified as: 

1. Topologically nontrivial thick branes: These solutions are based on the assumption of existence of some topo- 
logical defects in spacetime. Solutions in this class are composed of thick branes made of cither a single scalar 
field or non- interacting multi-scalar fields (and so on). 

2. Topologically trivial thick branes: These solutions can exist for interacting scalar fields but they are absent in 
the case of one scalar field. 

Several representative approaches to formulate cosmological equations on thick branes will be reviewed (section IVT). 
Finally, the regular S-brane solutions will also be explained (section IVl)) . 

As a final remark, we will not discuss applications of thick brane models to high-energy physics and the problem 
of localization of field with various spins. In this review, we are going to focus on their geometrical and topological 
properties only. 

II. TOPOLOGICALLY NONTRIVIAL THICK BRANES 
A. Definition of topologically nontrivial solutions 

In this review, we use the classification of thick brane solutions given in section ITD] According to this classification, 
one type of solution is topologically nontrivial. In this section we give the definition of topologically nontrivial 
solutions. This classification is based on definitions of kink-like and monopole-like solutions. 

Kink- like, topologically nontrivial solutions are considered in two dimensions (one space and one time dimensions) . 
Such solutions come up in considerations of problems with one scalar field (j) with a potential energy having two or 
more degenerate vacuum states. One example is the well-known Mexican hat potential (52j 

where A and are coupling and mass parameters respectively. In this case, there exist localized solutions which go 
to 4>±oo = ±to/-\/A (the so-called kink) asymptotically. The solution with an asymptotic behavior 0±oo = ^m/y/X is 
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known as the anti-kink. In the papers reviewed below the asymptotic values of the scalar field (j)±oo can be shifted by 
(/)o , i.e. (j)±oo = ±rn/VX + cj)o. 

We now consider the case when the dimensionality of space is n > 1 , and when the number of scalar fields is m > 1 . 
We take the following mapping of a sphere Sn-i on a sphere Sm-i- sphere Sn-i is an infinitely remote sphere in a 
space of extra dimensions, i.e. it is a sphere at ?' — > oo in the metric 

dll, = bir)idr^+r^dnl_,). (18) 

Here dl^j^ is the metric in the extra dimensions; the metric b(r)r^dfl^_i is the metric on a sphere Sn-i- The sphere 
Sm-i is a sphere in the abstract space of the scalar fields 0", a — 1, 2, • • • m. 
Using this set up the mapping is defined by 



n = > 

r 



A picture of this mapping is given schematically in Fig. [51 This mapping maps the point 
the scalar fields to the vector n = - G Sn-i in a space of the extra dimensions. 



(19) 

G Sm~i in the space of 




FIG. 2: Map n : 



for a case n ; 



From the topological point of view, we have the following map 

TTn-l {Sm-l) = Sn-1 ^ Sm-1, (20) 

where TTn-i (Sm-i) is a homotopy group of the map Sn-i — > Sm-i- It is known from topology that all nonsingular 
maps from one sphere Sn-i into another one Sm-i can be divided into homotopic sectors. Maps inside of one sector 
can pass continuously each into other, but maps of two different sectors cannot pass continuously each into other. 
For example [HSi 

^„(5„) = Z, (21) 
T^niSm) =0 by n<m, (22) 
^„(5i) =0 by n>l, (23) 

where Z is a group of integer numbers, i.e. the maps 5,1 — > Sn are divided in to a discrete set of homotopic classes 
which are characterized by integers. In ([22)) and ([23)) . zero on the right hand side means that the group is trivial, i.e. 
in this case all maps can be deformed each into other. 

We will call a thick brane solution topologically nontrivial if the corresponding homotopy group TTji—i {Sm—l^ 
nontrivial. Note that the case n = 1 and m = 1 is topologically nontrivial. 
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B. Thick brane solutions with a single scalar field 

111 this subsection, as the simplest example of a topologically nontrivial solution, we review the five-dimensional 
thick brane solutions with a single scalar field. 

1. General properties 

First, we review the general properties of thick brane solutions of a single scalar field which have the four-dimensional 
Poincare symmetry, following Ref. (53 . Issl ] . We consider a five-dimensional Einstein-scalar theory 

S=^Jd''xV^[R-{dcl,f~2V{(t>)]. (24) 

In general, the metric of a static five-dimensional spacetime with the four-dimcnsional Poincare symmetry can be 
written as 

ds2 ^ e2F(z)^^^^^;.^^;. _^ e«^(^)dz2. (25) 

For simplicity, it is assumed that the bulk geometry is Z2-symmetric across the center of the brane, hence F(—z) = 
F{z), 4>{—z) = —(j){z), and the imposed boundary conditions are 

F(0) = i^'(O) = 0, <?i(0)=0. (26) 

In Ref. [53, , the properties of thick brane solutions were investigated and it was shown that the only possible 
asymptotic geometry of the thick brane spacetime is asymptotically AdS, and is only possible if the potential V{(j)) 
has an alternating sign and satisfies the fine-tuning condition 

V{oo) = 0, V{z) / y/^V{(P(z))dz = / e^-^(^V(0(z))dz. (27) 
J a Jo 

The zero thickness limit is well-defined and the solution reduces to the Randall- Sundrum model if the asymptotic 
value of V{(j>) approaches a constant, which is the bulk cosmological constant. Note that if the brane has the four- 
dimensional de Sitter or anti-de Sitter symmetry, the fine-tuning condition (P7| may not be satisfied. 

2. Single scalar thick brane solutions 

In this subsection, we review the various thick brane solutions whose geometry can be Minkowski, de Sitter or anti- 
de Sitter, although in terms of the cosmological applications we focus on the first two cases. Many of previously 
known solutions arc discussed in the context of the Z2-symmctric bulk spacetime but non-Z2 symmetric solutions 
have also been reported. The thick brane solutions of more than two scalar fields or of other fields will be reported, 
separately. 

a. Set-up We focus on the Einstein-scalar theory given in Eq. (P^ . Most of solutions reviewed here have the 
following form of the metric 

ds^ = dy^ + a(y)27^,c^a;''dx^ = ^(y) (28) 

with the four-dimensional Ricci curvature tensor being -R[7;.ii/] = 12K. The constant K — 0, 1, —1 represents the four- 
dimensional Minkowski, de Sitter and AdS sections, respectively. It is also assumed that the scalar field depends only 
on the extra space coordinate y: (f> = (t>{y). If the Z2-symmetry is imposed, namely a{y) = a(— y) and (t>{y) = —(j){—y) 
(the center of the brane is placed at y = 0). the boundary conditions for metric and scalar field arc given by 

a'(0) = 0(0) = 0, (29) 

where "prime" denotes a derivative with respect to y. The warp factor is normalized so that a(0) = 1. At the infinity 
y — >■ ±00, the scalar field approaches some constant value 4> = ±0oo, where for a kink > and for an anti-kink 
0OO < 0. By minimizing the action Eq. p4)). the Einstein and scalar equations are explicitly given by 
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Sometimes, it is useful to reduce the above Einstein-scalar equations into first order equations. This can be achieved 
by introducing the auxiliary superpotential W{<j)) jSGj - lsTj . in terms of which the scalar field potential is given as 



V = -6M^(^)2 + 



9 /dW 



7:= \ 1 



K 



Using this, the above equations become 



3 dW a' 

J dtp a 



(31) 



(32) 



For a given superpotential W the scalar field profile is obtained by integrating the first equation. Then, the metric 
function is determined by integrating the second equation. 

More directly, for a given metric function, one can determine the scalar field profile and potential, using the 
procedure of [5l|. But sometimes, it is impossible to give the potential as an analytical function of the scalar field, 
rather than as a function of the bulk coordinate. Another method [111, [G^I is that first one gives the profile of the 
scalar field, and then the metric function and the potential are determined through the field equations. 

b. Minkowski brane solutions Here the simplest cases are solutions which contain the four-dimensional section 
which has the Poincarc symmetry K = 0. The solutions which are obtained from the simple superpotential 



W = csin(5(/)) 



(33) 



were discussed in [57H59l . l63l - l65l | (See also the textbook [50[). The corresponding scalar potential has the sine-Gordon 
form 

3 



V{(j>) ^ -c^ [3&2 cos^ (60) -AsiT? 



The solution is given by 



cosh(c6^y) 



1/36^ 



■ arctan 



tanh 



ficb'^y 



In the asymptotic limit y — !• 00, 0(±oo) = ±7r/(2&) the potential approaches the negative values V{±oo) 
The bulk is asymptotically AdS. This class of solutions also contains the solution derived in [5l| . 



(34) 
(35) 



a{y) 
0(y) 



1 



^-■2n(y+yo) _|_ g2n{y+yo) 

±\ - arctan(e2"''), 
V n 



l/2n 



whose potential is given by 



cos 




V{(p) = -6cos2 




3n sin 



6n 







(36) 



(37) 



The dimensionless parameter n represents the brane thickness. The solution is stable against perturbations. 
In Ref. (66j , another type of thick brane solution was discussed, which was obtained from the superpotential 



av 



1 



3^2 



leading to the potential 



A 

4 



271/2 



V 

"3 



(38) 



(39) 



where A :~ 2a^/v'^ is scalar self-coupling constant and the gravitational scale M is shown explicitly. In the absence 
of gravity A/ — > 00, the theory reduces to the Xcfi'^ scalar field model. The scalar and metric functions are given by 



(j>{y) = i;tanh(ay), 
g 

In a{y) = Cq - —v^ In [cosh(ay)] 



54 



tanh^(au). 



(40) 
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where Co is an integration constant which represents the overall scaling of our four-dimensional world (See also the 
textbook [131 )• It is trivial to check that the geometry reduces to the five-dimensional Minkowski geometry in the 
decoupling limit AI ^ oo. In Ref. the potential 



was considered. Here Vq > 0. This potential gave rise to thick Minkowski brane solution with the asymptotically AdS 
bulk. 

c. de Sitter brane solutions Here, we review some known thick dc Sitter brane solutions with K — +1. 
In Ref. [HH, the stable de Sitter brane version of the solution Eq. (|36ll was discussed, whose metric function is 
given by 



1 



-1 l/2n 



sinh ^"(j/ ?/o) + sinh^"(7y + yo) 



(41) 



But there is no analytic form for (/)(?/) and V{4)) |5l| . 
In Ref. [13, HH , for the scalar field potential 



9i?2(/?2 _ 1) 15i72(^2 ^ 



2/32 



2/32 



cos 



(42) 



the stable solutions 



-sn(iJy,*/3-i), 



(f>{y) ~ VG arctan 
were derived, where the elliptic functions are defined by 

sn^"'"(z, k) 



cn(77y, i/3 



_l3dn{Hy,il3- 



dx 



^{l-x^){l - fc2a;2)' 



(43) 



(44) 



cn(z) = (1 - sn2(w, k)Y^'^, and dn(z) = (1 - k'^sn'^{u, k))^/"^. 

Other interesting solutions were given for an axionic scalar field potential [gsI-ItS 



F(0) = Focos2(i--) -f 



(45) 



where < ct < 1. The Z2-symmetric solution [68| is given by 



a{z) 



cosh 



Hz 
a 



= (f>o arcsin 



tanh 



Hz 
a 



(46) 



where dy = a(z)dz and the brane expansion rate H is related to the thickness and potential parameter through 

2aVQ 



H' 



3 l + 3cr 



(47) 



The physical thickness of the brane is given by a/H. In the limit tr ^ 0, the solution smoothly approaches a thin de 
Sitter brane embedded in a five-dimensional Minkowski bulk. 

Extensions of Eq. pS)) to the non-Z2 symmetric bulk were considered in Ref. [7ll-[73j. Other thick de Sitter brane 
solutions with more complicated metric functions were reported in Ref. [6^ . 



C. Thick brane solutions in generalized scalar field theories 

In this subsection, we review the thick brane solutions in more general theories of a single scalar field, i.e., a scalar 
field with a non-standard kinetic term (t^ - It?! , a scalar field non-minimally coupled to gravity [78| , a phantom scalar 
field [Zi HO], and a tachyon field [U. 
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1. Scalar field with a non-standard kinetic term 

The paper [tJ discussed the case of a scalar field with non-standard kinetic term. Such a model was initially 
considered as the kinetic generalization of the single field infiation models [12, HI] and extensively studied for modeling 
inflation and dark energy. Here it is applied to a fivc-dimcnsional space-time to find new thick brane solutions. The 
starting point is the action 

S = j <fx^{^ - ii? + £(^, X)) , (48) 

where gab is the five-dimensional metric with the signature (+,—,—,—,—), and the five-dimensional gravitational 
constant is chosen G^^^ — l/(47r). The indices M, N, - ■ ■ and ii,i>, ■ ■ ■ run over the five- and four-dimensional space- 
time, respectively. We now define X {l/2)g^^^V m4'^ n4>- Then, the equations of motion are given by 

G^'^'dMdNd^ + 2XCx4> - = , (49) 

where G*^^ = Cxg^^^ + -Cjcx V*^0V^0. The Minkowski thick branc solutions whose metric is given by 

ds^ = e^'^'^y^^^^dx^'dx" - dy^, (50) 

are considered. The equations of motion are 



(' 



Cx + 2XCxxM" - hxCx^ -CA^ -ACxcji'A' 



^" = \xCx. 

A')^ = l^[c-2XCx), (51) 

where the prime is derivative with respect to y. For convenience, the (superpotcntial) function W{(l)), such that 
Cx<j)' = Wtf,; is introduced. The scalar field Lagrangian is assumed to be 

C = F{X) - V{4>) , (52) 

where F can be an arbitrary function of X. Note that F{X) = X for the standard canonical scalar field. The 
equations of motion now reduce to a set of the first order equations that 

F'(X)0' = ]^W^, 

F -2F'{X)X -V{(t>) ^\w^. (53) 

In Ref. [zll, two specific examples of the choice of the function F were considered: (I) Fj = X + and (II) 

Fn = -X\ 

a. Case (I) In case (I), Fj = X + a\X\X , it is difficult to derive solutions for general values of a. Thus, Ref. 
[t^ considered the case of small a and discussed perturbations around the case of a canonical scalar field. Then, the 
behaviors of scalar and metric functions are given by 

<t>{y)^My)^l<l^'^{y)w{My)). A(y)^ A^{y) + ^w{Uy))\ (54) 

where <?!'o(y) is the solution in the case of a = (the canonical case). For the specific choice that 1^(0) = 3asin(60), 
one finds the scalar potential 

V = ^a^fo^ cos2(60) - 3a2 sin^{bc^) - ^a%'^ cos\b4>) + 0{a^) . (55) 
8 64 

The scalar and metric functions were found as 

4>{y) = T arcsin tanh ( -ab^y ) \ 4- + 0{a ), 



2 , / , /3 



9aa'^b 


sinh 


4 


cosh^ 






tanh^ 


4 





For larger value of a, the peak of the warp factor A is spread out and the energy density at the center of the brane 
decreases. 
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b. Case (II) In case (II), Fu = —X'^, the choice of the superpotential, 

Wicj)) = 9a^b^ sin(6(/)) (2 + cos^ {b(j})j , (57) 

leads to 

y(0) ^ '^a^'^ cos\b^) - 27a**64 sin2 (60) (2 + cos2 (&(/))) ^ (58) 

An exact solution was found 

(kiy) = ^ arcsin (tanh(3a62y)) , A{y) = -y tanh^ (Safety) - ^a^ In cosh(3a6^y)^ . (59) 

The effective brane thickness is of order (a6^)^^, which can become larger for smaller a and 6. 

The work of [t^, [tI] considered the model composed of the non-linear kinetic term F{X) = X\X\ and the self- 
interacting potential V{4>) = A^((/)^ — (/)g)^, where cf) = ±(/)o corresponds to the global minima and A is the self-coupling 
of the scalar field. A solitonic solution coupled to gravity, which can be interpreted as a thick brane, was obtained 
by employing a series expansion method and a numerical one in [76j . while an exact solution without gravity was 
originally derived in (TSj . The solution is stable against the linear perturbations. In addition, this model has an 
appealing property that the propagation of perturbations of the scalar outside the thick brane is completely suppressed 
and the scalar field is automatically restricted to the brane. 

The recent work of [tJ discussed the de Sitter or AdS thick brane solutions in more general class of models that 

C = Fix, 0) - V{^) . (60) 

The line elements of the de Sitter and AdS brane solutions are assumed to be 

and 

ds^ = e'^^y^e-^^^'idf - dxl - dxi,) + dxl) - dy\ (62) 

respectively. The choice of F{X, (p) — X(t>™ , where m = 0, 1, 2, • • • is considered. The solutions are given in assuming 
that A = \n (cos(6y)). Note that in this ansatz there are naked singularities ai y — ±y* := ±7r/(26). 
The de Sitter brane solutions are obtained as follows: For m ^ 2n where n = 0, 1, 2, • • • , 

0(y) = [{n + l)/3arcsinh(tan(6y))] , (63) 



where /3 = i^(3/2)(fe2 - A), for the potential 



9, 



F(<^) =362-- (62 -A) ^Qgj^: 



(n + l)/3 



(64) 



Note that for an odd n the solution of cf>{y) is invalid for a negative y, while for an odd n the solution of (f>{y) is valid 
everywhere in —y* < y < y* . For m = 2n -I- 1 where n = 0, 1, 2, • • • , the de Sitter brane solutions are given by 

= [in + i)/3arcsinh(tan(6y))]'/^'"+'\ (65) 



where ^ = l^/{S/2){b'^ - A), for the potential 

-2^(2n+l)/2 



y(</)) =362-^(62-A) cosh^ 



(2n-f 1)/? 



(66) 



The solution is invalid for a negative 0. The de Sitter brane solutions exist only for < A < 6^. The AdS brane 
solutions are obtained by reversing the sign of A. Then, A must take a value in the range of — 6^ < A < 0. 

The case of m = recovers the solutions with naked singularities in the bulk obtained in the case of a canonical 
scalar field [H . 
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2. Non-minimally coupled scalar field 
Ref. [tsI considered thick brane solutions of a scalar field non-mininially coupled to the scalar curvature: 



S = / (^'x^/-G 



f{4>)R-\[d4)'' -Vicjy)], (67) 



where /(</)) denotes the direct coupling of the scalar field to gravity and the five-dimensional metric Gab has the 
signature (— , +, +, +, +). The above action is conformally related to the Einstein frame action with the Ricci scalar 
term {2M^R) via the conformal transformation Gab — > Gabf{4')/{'^^^): where AI is the fivc-dimensional gravitational 
scale. The metric ansatz was taken to be 

ds^ = e^'^rj^^dx^'dx" + dy^. (68) 

The coupling function /(</>) was chosen to be /((/>) = 2A'P — \£,(f)^ ^ where the cases of ^ = and ^ — 3/16 represent 
the minimal and conformal couplings, respectively. The Einstein equations are given by 

y (0) = - ^ {2M^ - ^4>^) (2i2 + i) + + ^ 

^ - |</,2) A - + i'cjycp + . (69) 



In the examples below, the scalar field potential is obtained once the scalar and metric profiles are found. 
In the simplest case, ^ = 0, the solution was given by 

(/.(y) = 0otanh(a2/), = (cosh(ay))"^e-^ (70) 

where 7 = (/)q/(9M'^). For a non-zero coupling constant ^ 7^ 0, by choosing 



h(ay))^^ (71) 



the solution 4i{y) = 0o tanh(a2/) was obtained with 



The above solution exists for < ^ < 1/6. The scalar curvature of this solution approaches R — 20a^(^ ^ — 6) 
for ?/ — > 00, namely the bulk is asymptotically AdS. Eq. ([7T|) is also the metric function of the solution whose scalar 
field profile is given by (/)(y) = 0(O)( cosh(ay))"\ with 0^(0) = 12M^{^-^ - 6)/(3 - 160- 
Another solution with 4>{y) — 0otanh(aj/), where 4>o = 2y/S,^^M^, and 

A{y) = -4 In ( cosh(ay)) + ^ (8 - T tanh^ (ay )^^pfQ ({l, 1, 7/6}, {3/2, 2}, tanh2(a2/)) , (73) 

where Fppq represents the generalized hypcrgcomctric function, was also found. Note that in all the above nontrivial 
solutions there is no smooth limit to the minimally coupled case ^ = 0. 



3. Phantom scalar field 



Ref. [79| considered a model of a bulk phantom scalar field: 

1 



S 



d^Xy/^ 



24 



R-2A 



(74) 



where gab is the five dimensional metric with signature (— , +, +, +, +) and V{(j)) is the potential of the scalar field (/). 
Note that sign of the kinetic term is now negative, i.c, phantom. The thick brane solutions for the phantom scalar 
field with the ansatz 



ds2 ^ ^ e''^f^y^r|^,^dx^'dx' 



(75) 
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and (j) = 4){y) arc considered. The equations of motion are given by 



A 



6 ' 

(76) 



/ del) 

where a :— k|/3. An expUcit solution was obtained for the sine-Gordon potential 

y(0)=B(l + cos(^)) (77) 




(78) 



where ki := (l/A^) ^|A|/6. Here the constant B has to be chosen as i? = (|A|/6a^) ya — 1/(4A^) j . The solution has 

a growing warp factor for the larger value of \y\, in contrast to the case of the normal kinetic term. 

In [sot , by using a massless scalar phantom/ghost field, an analytic time-dependent warped solution forming a 
standing wave in the bulk was found. The nodes of the standing wave are interpreted as the different four-dimensional 
Minkowski vacua, called islands in Ref. [lOl , having different physical parameters such as gravitational and cosmolog- 
ical constants. This model possesses the following characteristic features: 1) The ordinary four-dimensional matter 
can be localized on the tensionless branes, i.e. on the nodes. 2) It is possible to realize a new gravitational localiza- 
tion mechanism when matter fields can reside only on Minkowski islands, where the background space-time does not 
oscillate. 



4-. Tachyon field 



by 



As a final example of a generalized scalar field theory, consider a model of a tachyon field [8 1| whose action is given 



(79) 



where A5 and V(T) are the bulk cosmological constant and the tachyon potential, respectively, g is the determinant 
of the five-dimensional metric gab whose signature is (— , +, +, +, +). It is assumed that the metric and tachyon field 
take the following form 



ds^ = dy"^ + e^f^y^rj^^dx^dx" 



(80) 



and T = T{y). 

A thick brane solution is obtained for non-zero cosmological constant A5 ^ 0. The potential is explicitly given as 
a function of y, as 

V{y) = \ , 2 ^ \/(A5 + 66^) sinh^(b7/) + A5 + 6(6^ - m)J{k^ + 66^) sinh^(b7/) + (A5 - 6g^) , (81) 
K5 cosh [byj ^ ^ 

where h is an arbitrary constant. The warp and tachyon functions are given by 

f{y) = lncosh(6y), 

' Ag + 662 



ny) 



i_ Sjb^ + H^) 
by A5-6iI2 



Elliptic^ 



iby, 



A - 6i?2 



(82) 



The warp factor grows for large y. 
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D. Other types of thick branes 

1. Multi-scalar 

The brane models with two interac ting: scalar fields considered in section UlI Al refer to the so-called non-topological 
solutions in the terminology of Ref. [S^]- At the same time, there exist other types of solutions - topological ones. 
These solutions describe our 4-dimensional Universe as a global topological defect in a multidimensional spacetime. 
Such a topological defect can be created by a set of scalar fields describing some "hedgehog" configuration (85l - [87| . 

The main idea of the research in this direction is illustrated by the example of the paper [85| where the {p— l)-brane 
solutions (the number p refers to the coordinates on the brane) are considered. It is assumed that there is a global 
defect in n extra dimensions which is described by a multiplct of n scalar fields, 0°, whose Lagrangian is ^ 

L^-\dArd^r~V{cl>), (83) 

where the potential V{(j)) has its minimum on the n-spherc — rf. As an example, the potential 

v{<t^) = -^{rr-i^f (84) 

was used. Using the following ansatz for the metric: 

ds^ = ^(O'^e' + C'd^^'-i + B{ifg^,dx''dx\ (85) 

where dVl^^ stands for the metric on the unit m-sphere, and the spherical coordinates in the extra dimensions are 
defined by the usual relations, ~ {^cos^i, ^ sin 6*1 cos 02, •■•}■ 

Using a different metric ansatz, a few different types of solutions were found. The first typ e, without a cosmological 
constant, exists for all n > 3, and is very similar to the global monopole solutions from [88|. The brane worldshcet is 
flat, and there is a solid angle deficit in the extra dimensions. 

For a positive cosmological constant, A > 0, the solutions describe spherical branes in an inflating multidimensional 
universe. In the limit 77 — > 0, when the gravitational effect of the defect can be neglected, the universe can be 
pictured as an expanding (p + rt — l)-dimensional sphere with a brane wrapping around it in the form of a sphere 
of lower dimensionality (p — 1). Another class of solutions has curvature singularities even in the absence of a defect 
(77 = 0), and these solutions were considered as unphysical. Finally, there are solutions having the geometry of a 
ij) + l)-dimensional de Sitter space, with the remaining [n — 1) dimensions having the geometry of a cylinder. 

For a negative cosmological constant, A < 0, three classes of solutions were found. The first two are essentially 
analytic continuations of the positive-A solutions. The third class is similar to the Randall-Sundrum [n = 1) and the 
Gregory (n = 2) solutions, and exhibit an exponential warp factor. 

In [83| a model of a global topological defect with the potential was considered. The authors were looking for 
solutions decaying far from the brane as B{£) ~ exp(— c <^). It was shown that there exist solutions with the spherical 
radius R const at large ^, i.e., the extra dimensions form an ?7-dimensional cylinder x S*""^. Such solutions 
exist only for n > 3. 

The branes considered above refer to Minkowski branes with some non-zero bulk cosmological constant. In the paper 
[sot a model with some extra 4-dimensional A-term is considered (i.e. a de Sitter brane). This work investigated global 
monopole solutions in the Goldstone model with n scalar fields with the Lagrangian ([SS]). Within the framework of 
this model, it was shown that new types of solutions exist (mirror symmetric solutions) . These solutions have periodic 
matter and metric functions, and require a fine tuning of two cosmological constants. 

2. Vortex 

If one replaces the domain wall with some more complicated topological defect, such as string- vortices or monopoles 
(see above), new possibilities for gravitational localization to the brane appear. Braneworld models using a 6- 
dimensional vortex were proposed in Ref. [l^l- In this model gravity was localized on an Abrikosov-Nielsen-Olsen 
vortex in the context of the Abelian Higgs model in 6 dimensions. In this model, the conservation of a topological 



^ In the original paper (85|| , the sign of the kinetic term is different but it is wrong. 
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charge guarantees the stabiUty of the configuration. In principle the presence of gravity can change the situation, but 
one can hope that "gravitating" topological defects can provide potential candidates for stable configurations. 
For this purpose, the total action of the gravitating Abelian Higgs model in six dimensions can be chosen as 

S ~ Sh (86) 

where S'brano is the gauge-Higgs action and S'grav is the six-dimensional generalization of the Einstein-Hilbert action. 
More specifically. 



^branc 



d^xT^/Ibrane, ^brane = i(PA0)*I?^0 " -FabF'''' - ^ (0> - 1^')' , (87) 



where T) a = a — i&Aa is the gauge covariant derivative, while V/i is the generally covariant derivative, and the bulk 
action is 

5grav = - / d^'x^TG (^+A 



I' 



2x 



where A is a bulk cosmological constant, x = SttGq = Stt/Mq and A/g denotes the six-dimensional Planck mass. In 
Eq. ([87)1 . V is the vacuum expectation value of the Higgs field which determines the masses of the Higgs boson and 
of the gauge boson 



rriH ~ v2Xv, my = ev. (89) 

Choosing the six-dimensional metric as 

ds^ = GAsdx^dx^ ^ A'P{p)g^,^dx^dx'' - dp^ - L{pfde^, (90) 

where p and 9 are, respectively, the bulk radius and the bulk angle, 5^1/ is the four-dimensional metric and M(p), 
L{p) are the warp factors. The Abrikosov-Nielscn-Olscn ansatz for the gauge-Higgs system reads: 

(/.(p,^?) = ^;/(p)e^"^ 

Ae{p,e)^-[n- P{p)\, (91) 
e 

where n is the winding number. 

Choosing the corresponding regular boundary conditions, one can find that outside the core of the string all source 
terms vanish and a general solution to the equations can be easily found for the case A < (solutions for the case of 
a positive cosmological constant in the bulk were studied in 91[) as follows 

m{x) = -c — , (92) 



where e is an integration constant; c = yJ—ji/lQ > 0; /i is some combination of the 6-dimensional Einstein gravitational 
constant, A-term and the mass of the Higgs; ■m{x) and x are the rescaled variables M[p) and p, respectively. Once 
m{x) is known, the function ^{x) (the rescaled L{p)) can also be determined. 

To clarify the question about gravitational localization on the brane, one can demand the finiteness of the four- 
dimensional Planck mass 

Ml = / dxM^{x)C{x) < 00 . (93) 

mjj J 

It was shown in that ife>Oore<— 1, then for a; ^ 00 one has M ^ C ^ ef^^ so that the integral in ([M)) 
diverges and gravity cannot be localized. 
If e = the solution is simply 

m{x) = l(x) = -c, (94) 
and the warp factors will be exponentially decreasing as a function of the bulk radius: 

M(x) = Moe-'=% 

Lix) = roe"'^'^. (95) 
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94 )) - ([95)) leads to the localization of gravity and to a smooth AdS geometry far from the string 
where = -r^- log 1 / 1 e 



The solution of Eqs 
core. 

If —1 < e < 0, A/(a;o) = where xq ~ -g^logl/lel. The geometry is singular at x — xq thus x < xq should be 
required. In spite of the fact that C diverges at x = xq, the integral (j93p defining the four-dimensional Planck mass 
is finite. So, these geometries can be potentially used for the localization of gravity, provided the sing ularity at xq is 
resolved in some way, e.g. by string theory. These types of singular solutions were discussed in [9l[ for the case of 
positive cosmological constant. 

Finally, if the bulk cosmological constant is zero, the solutions have a power-law behavior, namely 



M{x) 



Cix) 



with 



(96) 



(97) 



where d = 4 is the number of dimensions of the metric g^jy. These solutions belong to the Kasner class. The Kasner 
conditions leave open only two possibilities: cither S = 1 and 7 = or 7 = 2/5 and S = —3/5. None of them lead to 
the localization of gravity. 

Summarizing the results of the paper [g^l, it was shown that the localization of gravity is possible on a "thick" 
string, and a fine-tuning condition was found which led to a set of physically interesting solutions. Since the described 
geometries are regular, gravity can be described in classical terms both in the bulk and on the vortex. 

Similar considerations allowed the authors of [9^ to conclude that within the framework of an anomaly-free Abelian 
Higgs model coupled to gravity in a 6-dimcnsional spacctime, it is possible to construct an effective Z? = 4 electro- 
dynamics of charged particles interacting with photons and gravitons. In this work both gauge field and gravity are 
localized near the core of an Abrikosov-Nielsen-Olsen vortex configuration. 

Further considerations of questions concerning localization of gauge fields on a gravitating vortex in a 6-dimensional 
spacetime can be found in the papers [sl, [13] . Note that most of the calculations are performed numerically or 
asymptotically, and there are no general analytic solutions for 6-dimensional vortices. However, in order to address 
some questions (for example, stability analysis) the availability of exact solutions is desirable. In Ref. [9^ some exact 
solutions for 6-dimensional vortices were found. In particular, using the special boundary conditions, the several kink 
solutions were found. The first solution was 



/ = /o arctan ^sinh , fo = 2V5, 
with the corresponding metric "warp" factor 

S>0, 13 >0, 



M{r) = cosh-^' ( ^ 



(98) 



(99) 



where 6 could be seen as the wall thickness and /3 as a parameter depending on the cosmological constant. The 
effective potential for this case was 



n/) = ^[(l + 5<5)cos2 (///„) -55]. 

The second and third solutions were, respectively 

f = fo tanh (ar) , and 
f ^ fa arctan (ar) , 

with the metric "warp" factors, respectively 

cosh~-''° /^(ar). 



M(r) 
M{r) 



exp 



exp 



^/o sech(ar) 



and 



-/qT arctan(ar) 



(100) 

(101) 
(102) 

(103) 
(104) 



and with a polynomial potential identical to the Mexican hat potential, or a potential containing polynomial and 
non-polynomial terms 



V{f) 



a 

y 



{f'-foY'-lf {f-ifi) 



, and 



(105) 



20 



Using these exact solutions, it was shown that the zero mode of the linearized gravity spectrum is localized on the 
3-brane and the massive modes are not bounded at all. 

3. Weyl gravity model 

In Refs. (96l - [98j . thick brane solutions in a pure geometric Weyl 5D spacetime, which constitutes a non-Riemannian 
generalization of Kaluza-Klein gravity, are considered. A Weyl geometry is an affine manifold specified by (gMNj^M) 
with M,N ~ 0, 1,2,3,5, where qmn is the metric tensor, and ujm is a "gauge" vector involved in the definitions of 
the affine connections of the manifold. 

Besides the scalar curvature R, the action used in this approach also contains terms describing a Weyl scalar uj 

I "^^r^ ei^fi? + iaViof + W{u)l (107) 



IGTrGn 



rW 



RmN — ^MN,A ~ ^\m.N + rf/TyrpQ — rf/gF^p, (108) 



where Mg is a Weyl manifold specified by the pair {gMN,i^)- The Weylian Ricci tensor reads 

MJV^ PQ ~ ^ MQ^ NP^ 

where 

Ta/jv = {mn} - \ (^,M<5^ + ^m5Z - gMNio'^) (109) 

are the affine connections on , arc the Christoffel symbols and AI,N ~ 0,1,2,3,5; the constant is an 

arbitrary coupling parameter, and U{uj) is a self-interaction potential for the scalar field u. 
The 5-dimensional metric describing a thick brane is assumed to have the form 

dsl = e^'^^y'>r]mndx"'dx'^ + dy^, (110) 

where e^^^^^ is the warp factor depending on the extra coordinate y, and m,n = 0, 1, 2,3, rjmn is the 4D Minkowski 
metric. 

Further, the conformal transformation qmn ~ e'^ Sain, mapping the Weylian action (|107p into the Riemannian one 
is carried out via 



^5 

M, 



167rG5 



R + iiiVujf +&U{u) , (111) 



where £,=£, — I, U{uj) = e '^U (w) and all hatted magnitudes and operators are defined in the Riemann frame. Now 
the metric (jllOp takes the form 

dsl = e^''(^)ry„™dx"dx" + e^'^^^dy^, (112) 
where 2a = 2 A + w. If one introduces new functions X = uj' and Y = 2 A' [o^ then the Einstein equations read 

X' + 2YX + -X^ = -—e'^, (113) 
2 (, duj 

Y' + 2Y^ + -XY= (--—+4U]e^. (114) 
2 \ ^ / 

As pointed out in [9^, this system of equations can be easily solved if one uses the condition X = kY, where k is an 
arbitrary constant parameter, but excluding the value k = 1. 

In this case both field equations in (|113p (jll4p reduce to a single differential equation 

y' + i±^y2^ J^e(^H+i)", (115) 
2 1+fc ^ ' 

where A is some constant parameter. The authors of [96j first considered the problem for a case of Z2-symmetric 
manifolds. Here, the two following cases are considered: 

(i) ^ = —(1 + fc)/(4fc) and leaving k arbitrary except for k = —4/3, 

(ii) k = —4/3 and leaving ^ arbitrary except for ^ — —1/16. 
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The case (i) For this case, a solution was found for the simphfied (jllSp in the form 

uj[y) ~ 6fc ln[cosh(ay)], (116) 
[cosh(a2/)]^ (117) 



e 



2A(-y) _ 



where 



3k 



The case (ii) The authors of [971 . lOq considered another simphfied case with /c = —4/3 and ^ an arbitrary 
parameter. Then Eq. (|115p takes the form 

r' + 12AeP'^ =0 or lu" - WXeP"^ = 0, (118) 

where p = 1 + 16^. Note that the choice ^ = —1/16 leads to a constant potential and thus the case ^ 7^ —1/16 is 
assumed. For such a choice of the parameter, one can find the following solution for lu and A 



^ _2jn(^^^cosh[ci(y-C2)] !>, (119) 



^2A ^ }}LJ^cosh[c^{y^C2)]} . (120) 



Cl 

The following physically interesting problems were considered in detail: 

• A > 0, p < 0, Cl > 0. In this case the solution is a thick brane, and the fifth coordinate changes within the 
limits —00 < y < 00. 

• A > 0, p > 0, Cl = iqi. In this case the solution describes a compact manifold along the extra dimension with 

< (7i(y — C2) < TT and consequently this solution is not a thick brane. 

If one calculates the function uj'{y) for both the cases (i) and (ii), then one has 

Lu'iy) = abtanh{ay), (i) (121) 

iu'iy) = -^tanh[ci(y-C2)], (ii) (122) 
P 

respectively. As one can see, both solutions are kink-like. Thus these solutions are topologically nontrivial solutions. 

The de Sitter thick brane solution in Weyl gravity was investigated in Ref. [9^ . To show the analytic form of the 
solution, it is convenient to transform the fifth coordinate to the conformal coordinate z, defined by dz = e^^^'^^dy, 

dsl = e2^(^) {-dt" + e2/^*dx2 + dz'') , (123) 

where j3 is the expansion rate of the de Sitter spacetime. For the potential 

1 + 35,. 



U{uj) = ——l5'e-^[cos—\ , (124) 
2b \ ojQ/ 



where ujq = \/3S{l — 5) (0 < (5 < 1), the solutions are given by 

Lu^LUo arctan fsinh , e'^ = cosh-^* (^) e^^" ^'^'^'^ i^) . (125) 



After transforming back to the proper coordinate y, the geometry is asymmetric across the y = surface. But the 
solution still has a kink-like profile and connects two Minkowski vacua. 
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4- Block brane 



It is possible to borrow ideas from condensed matter physics and apply them to brane world models. Of interest 
is the situation which occurs in ferromagnetic systems, in which one has an Ising or Bloch type wall. An Ising wall 
is a simple interface with no internal structure whereas a Bloch wall is an interface which has a nontrivial internal 
structure. One can find thick brane solutions with internal structures like that of the Bloch wall. 

In the paper jlOO| . the thick brane solution with two scalar fields (</>, x), called a Bloch brane, was considered. The 
solution has an internal structure because of the dependence on x- For this case, the action is taken to have the form 



Sc = (i^xdy\ 



1 



1 



The thick brane metric is 



ds^ — gabdx'^dx^ 



(126) 



(127) 



where a, 6 = 0, 1, 2, 3, 4, and e is the warp factor, a; = y is the extra dimension coordinate. The potential is given 
by 



K(0,X) 



1 



-yv' 



where r is a constant. It is readily seen that this potential is unbounded below at x = const 

Vc((p,x) — > -oo. 
The corresponding equations describing two gravitating scalar fields are 

dV{cj>,x) 



4>"+AA'(t)' = 

x" + 4A'x' - 

A' = 

A'^ = 



dVi(l),x) 
dx 

-I (*"^ 



1, 



6 



+ x'^)-^v{<i>,x), 



where prime stands for derivative with respect to y. One can reduce these equations to the form 

1 dWr 



2 d(j) 
1 dWr 



X 
A' 

a. Bloch walls The solution of these equations is 

Hy) = tanh(2r2/), 

x{y) 



2 dx ' 

1 

3 



9r 



— 2 scch(2ry). 
[(1 - 3r)tanh^(2ry) - 2 lncosh(2r?/)] 



(128) 
(129) 

(130) 

(131) 
(132) 
(133) 
(134) 

(135) 
(136) 
(137) 

(138) 
(139) 
(140) 



The brane has a typical width as large as 1/r. In the case r < 1/2, the two-field solution represents a Bloch-type 
brane, where (0, x) represents an elliptic trajectory in the field space. For r > 1/2, the two-field solution is changed 
to a one-field solution, i.e., an Ising-type brane, connecting two minima {(t>,x) = (iljO) straight along the line x 
in the field space. Equation p38p shows that this solution is a topologically nontrivial one. 
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b. General method In Rcf . |10l| , the method to obtain more general (often numerical) Bloch brane solutions was 
given. The starting point was the general class of the potential 



W^e(0,x) -20 



(141) 



which includes the previous case, for the choice that a = 1, A = — 1 and fi ~ — r. From Eq. (|137p . dtp/dx = T^c,0/W^c,xi 
where Wc,^ := dWc/d(l), one obtains 



A(0^-a^) + MX^ 



(142) 



Introducing the new variable p 



one can replace Eq. (jl42p with an ordinary differential equation for p with 



respect to x, which can be analytically solved as 



p{x) ^4>^~a'= cox^'^ - Y^X' , ( A ^ 2/i) , 
p(x)=<^'-a'=x'(ln(x) + ci), (A = 2^^), 



(143) 



where cq and ci are integration constants. 
(l/2)W^c,x, 



Substituting Eq. (|143p into the differential equation for x, dx/dy 



dx 
dy 

dx 
dy 



±2fixja^ + coX^/^' 



X-2p 



X^ (A^2/i), 



^=±2/.x^a2 + x2(ln(x)+ci), (A = 2^). 



For the metric function A, 



dA 
dy 



dA dx 
dx dy 



dy 



-Wr 



(144) 



(145) 



which leads to dA/dx = ^(l/3)(VFc/M^c,x)- This can be solved as 



^(x) = "0 

^(x) 



2a2 1 A 



2a' 



Ql 



lux 



3fi - 



- 3m 2 
^2fi^ 
Aci 



18m 



-X 



9^ 

X 2 

X 

D/i 



(A ^ 2m) 
(inx-^), (A = 2m), 



(146) 



where ag and ai are integration constants which are chosen so that A{y = 0) = 0. 

To summarize, the general method to obtain solutions is as follows: First, for a given set of parameters A, a and 
M, one solves Eq. (|144p and determines the profile of x as a function of y. Then, from Eq. (|143p and (|146p . one 
can determine the profile of (j) and A as the function of y, through xiv)- I^i general, Eq. (|144p does not have exact 
solutions except for the limited cases discussed below. Note that the solutions discussed in Ref. jlOOl | correspond to 
Co = 0. It should be stressed that the model admits a particular class of solutions which cannot be obtained from the 
method given above. A set of solutions such that x = (for which Eq. (|142p does not make any sense) do not have 
any internal structure for the brane, although the s ystem admits a solution given by ^ = ±atanh(Aa?/). 

According to the method shown above, in Ref. |lOl| . several new analytic thick brane solutions were found for 
Co 7^ 0, which were called degenerate Bloch walls and critical Bloch walls. Here, by "degenerate" we mean that one 
cannot specify the property of the solutions only by fixing the potential parameters A, a and fi, but one also needs to 
fix the integration constant co. In fact, the solutions still have remarkable variety for different values of cq. 

c. Degenerate Bloch walls For the case A = m and cq < —2, the following solution was obtained 



x(y) = 



2a 



4cosh(2Ma2/) — cg 



V4 



a/cq — 4cosh(2May) — co 



and 



2a 



-I (4a=)/9 



■\/cq — 4cosh(2/iay) — co 



exp 



2a (cq — cq^Cq — icosh{2p.ay) — 4a) 
9(-\/c§ — Acosh{2p.ay) — cq)^ 



(147) 



(148) 
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where N is chosen so that A{y = 0) = 0. 

For the choice X = 4^ and Cq < 1/16, a similar solution was obtained 



and 



x{y) = - 



N 



2a 



\/\/l — 16co cosh(4a^y) + 1 



-2a 



sinh(4/iay) 



\/ ^J\ — 16co cosh(4a/^?/) + 1 



- 16a^/9 




f 4a2 


exp < 


r 9 



•\/Vl — 16co"cosh(4a/i?/) + 1 ' 

1 + 32a2co + VI - 16co cosh(4^ay) 
(Vl — 16co cosh(4^ay) + l)^ 



(149) 



(150) 



and again N was chosen so that A(y = 0) = 0. An interesting feature of these solutions is that, for some value of Cq, 
exhibits a double kink profile. This represents the formation of a double wall structure, extended along the spatial 
dimension. 

d. Critical Block walls For the case that X = fi and cq = —2, the solution 



X 



1 ± tanh(/ia?/) 



tanh(/iay) =F 1 



and 



„2A 



N 



(1 ± tanh(/iaj/)) 



2a^/9 



exp 



a 



(1 ± tanh(/ia?/))^ (1 i tanh(/^a?/))^ 

a 



was obtained. Similarly, for the case that cq — 1/16 and A ~ Afi, another solution was obtained of the form 

r- cosh(/ia'(/) ± sinh(/ia'(/) , a ,^ , ,„ .n 

X = -V2a , \ <^=Tr(lTtanh(2^ay)), 

\/cosli(2/xa2/) ^ 



and 



N 



cos\v{2^ay) 



16a^/9 



X exp ■ 



2a2 
~9~ 



tanh(2a/xy) [a^ tanh(2a/i?/) T (1 + 2a^)] | . 



(151) 



(152) 



(153) 



(154) 



The novelty in these cases is the fact that both, and x fields have a kink-like profile and the warp factor has the 
remarkable feature of having two Minkowski- type regions. 



III. TOPOLOGICALLY TRIVIAL THICK BRANES 



A. Thick branes of two strongly interacting scalar fields 



Here, we review the topologically trivial solutions. For example, one possibility is to consider two nonlinear grav- 
itating scalar fields which can create a 4-dime nsional br ane in a multidimensional spa cctim e. Such problems have 
been investigated for t he 5- dimcnsional case in jl02l Il03l | , for the 6-dimcnsional case in [l04| , and for the 7- and the 
8-dimensional cases in |l05| . From the physical point of view, the situation is as follows: an interaction potential of 
these fields has two local and two global minima. Thus there are two different vacuums. The multidimensional space 
is filled with these scalar fields which are located at the vacuum in which the scalar fields are at the local minimum. 
Therefore there is a defect in the form of a 4-dimensional brane on the background of this vacuum. 

The general approach to such solutions is as follows: the action of = 4 -t- ji dimensional gravity for all models 
can be written as Il06ll: 



S 



-R + Lr, 



(155) 



where M is the fundamental mass scale, and n is a number of extra dimensions. The Lagrangian, i„ 
interacting scalar fields Lp, x takes the form: 



for the two 



(156) 
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where the potential energy V{ip, x) is: 



^1/2 2n2 , ^2/ 2 



2\2 



y + ip'x' - Vo 



(157) 



4 "'4 

(This potential was used in jl07l | as an effective description of a condensate of gauge field in SU(3) Yang-Mills theory, 
i.e. the scalar fields were taken as effective fields which described a condensate of Yang-Mills fields.) Here the capital 
Latin indices run over A,B = 0,1,2,3,...,D and the small Greek indices a, f3 = 0,1,2,3 refer to four dimensions; 
Ai,A2 are the self-coupling constants, TOi,m2 arc the masses of the scalar fields ip and X: respectively; Vq is an 
arbitrary normalization constant which can be chosen based of physical motivations; e ~ +1 for usual scalar fields, 
and e = — 1 for phantom scalar fields. 

The use of two fields ensures the presence of two global minima of the potential (|157p at = 0, x = ± m2 
and two local minima at x = 0, ^ = ± mi for the values of the parameters Ai,A2 used in the above papers. 
The conditions for existence of the local minima are: Ai > 0,mi > A2m|/2, and for the global minima they are: 
A2 > 0,m2 > Aim^/2. Because of these minima there were solutions localized on the brane for the 5-, 6-, 7- and 
8-dimensional cases [l02l |- |l05| . In these different cases the solutions go asymptotically to one of the local minima. 

Variation of the action (|155p with respect to the Z?-dimensional metric tensor gj^g leads to Einstein's equations: 



1 



1 



(158) 



2 " KP' 

where i?^ and are the D-dimensional Ricci and the energy-momentum tensors, respectively. The corresponding 
scalar field equations can be obtained from (|155p by variation with respect to the field variables x- These equations 
are 



1 



d 



9 



d{(p,x) 



dx^ 



dV 



(159) 



1. General equations 



Using the generalized _D-dimensional metric |10f 

ds^ = (j)^{r)r]ap{x'')dx"dx'^ - X{r){dr^ + r^dnl^^) , 
where dfl^^^i is the solid angle for the {n — 1) sphere, one can rewrite the Einstein equations jl06j in the form 



(160) 



3 2 



'\2 



+ {n - 1) 



b' /A' . 2\ . A" lA' /A' 6\ . n-4A' /A' . 4 

T " 2T V a" " 



A r 



12^ + (n-l) 



4 A V A r 
2A 



]\Jn+2 a 

,0' fx' , 2\ n-2X fx' 4 
r) ^ 4 T V a" ^ r 
2A 



VA 



(161) 



(162) 



4 2 



0" T 



12^ + (-2) 



A r 



A 2 A V A 







n-5A' 










+ 4 A 


(f 


^9] 



(163) 



2A 



M 



n+2 



J 



where T^,T^ and Tg are (") , {^) and (g) components of the energy-momentum tensor, respectively. The scalar field 
equations become 



x" 



71-1 d)' n - 2 A' 

h 4— H 

r (j) 2 X 

n-1 6' n-2X' 



-4— 



2 A 



= X^[2x^ +Ai{^^ -ml)] 
X = Xx[2ip^+A2{x^-ml)] 



(164) 
(165) 
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The set of equations (|161|) - (jl65p is a set of nonlinear equat ions with solutions whose behavior very much depends on 



values of the parameters toi, m2, Ai, A2. As shown in Refs. jl02l |- |l05l |. by specifying some values of the self-coupling 
constants Ai, A2, the problem of searching for regular solutions to systems similar to the system given by (|161|) - (|165|) 
reduced to the evaluation of eigenvalues of the parameters mi,m2. Only for specific values of these parameters one 
was able to find regular solutions with finite energy. 



2. 5- dimensional problems 



The usual scalar field case with e = +1 was considered in [l02{ . and the phantom field case with e = — 1 was 
considered in |l03l |. These problems were solved by the numerical shooting method. The resul ts a re presented in 
Figs.|3]|4l These solutions describe a thick brane in a 5-dimensional bulk. As pointed out in Ref. [l03| . there are two 
distinctive features of the solutions for the phantom versus the usual scalar fields: a) they have different asymptotic 
behavior; b) the solutions with phantom scalar fields are st able, and the solutions with the usual scalar fields are not 
(at least for the parameters Ai, A2, (^(0), x(0) used in [l03|). 
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FIG. 3: The scalar fields profiles for the cases e = — 1 (solid 
lines) and e = +1 (dashed lines). 5D case. 



FIG. 4: The profiles of the functions % — 4>' l<i> for the cases 
e = — 1 (solid line) and e = +1 (dashed line). 5D case. 



3. 6, 7, 8-dimensional problems 

Thes e case s were considered in the papers [I04ll05ll . The results are presented in Figs.[S][ni All the results obtained 
in |l02l . Il03| for the 5-dimensional case and in jl04| for the 6-dimensional case, and also in [l05l | for the 7- and 
8-dimensional cases show that it is in principle possible to localize scalar fields with the potential (jl57p to a brane in 
any dimensions. This lends support to the idea that there exist similar regul ar so lutions is possible in a larger number 
(any number) of extra dimensions (a discussion of this question is given in jl05| ). 



4- Phantom thick de Sitter brane solutions in multidimensions 

Thick de Sitter brane solutions for phantom scalar fields were considered in |l08l ] for 5-, 6- and 7-dimensional cases. 
For this purpose the following metric was used 

ds"^ = a^{r)-iaf}{x'')dx°'dx'^ - \{r){dr^ + r^dnl_^) , (166) 

where dVL^_^ is the solid angle of the {n — 1) sphere. 7^^ is the metric of the four-dimensional de Sitter space whose 
Ricci curvature tensor is given by R^i,[l\ — '^H'^l^iv Then one has Einstein equations similar to equations (|16ip - (|163p 
with the following extra terms on the left hand side: —QH^X/cfP in the („) component, and —12H^X/(f)^ in the {^) 
and (g) components. The scalar field equations remain the same. The results of numerical calculations are presented 
in Figs. mfTUl 
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FIG. 5: The scalar fields profiles for the 6,7,8-dimensional FIG. 6: The profiles of the metric functions (j) and A for the 
problems. 6,7,8-dimensional problems. 
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FIG. 7: The scalar field configurations ^ and x are shown 
as functions of the dimensionless r. 



FIG. 8: The metric function A(r) is shown as a function of 
the dimensionless r. 



As one can sec from these figures, there exist five-, six- and seven-dimensional thick de Sitter brane solutions of two 
interacting {phantom) scalar fields. It was also shown in |l08j that asymptotically there exist anti-de Sitter spacetimes 
for the five-dimensional case and fiat spacetimes for the six- and seven-dimensional cases. For the five-dimensional 
case, it was shown that the solutions are stable. It is also quite reasonable to expect that the thick brane solutions in 
multidimensions are also stable since the spacetime and vacuum structures are essentially the same as in the case of 
five spacetime dimensions. 



IV. BRANES WITH UNUSUAL SOURCE FUNCTIONS 



Here we present the thick brane solutions created by special source functions [l09l . lllCll | . It was shown that by going 
from 5D to 6D it is possible to trap fields of all spin s (i.e. spins 0, ^, 1, 2) to the 4D spacetime using only gravity. 
The "trapping" provided by the 6D solution of |110l | is automatic since the proper distance from the brane at ?' = 
is finite and any field can only be distributed within some finite distance from the brane. 

The 6D action was chosen as 



S = 



—R + A + L„ 



(167) 



where M, i?, A and L„i sltc the fundamental scale, the scalar curvature, the cosmological constant and the matter 
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field Lagrangian, respectively. All these physical quantities refer to the 6-dimensional spacetime with the signature 

of (H ). The variation of the action ()167p with respect to the 6-dinicnsional metric tensor gAB leads to the 

Einstein equations 

Rab - ^QabR = (5abA + Tab) , (168) 

where Rab and Tab are the Ricci and the energy-momentum tensors, respectively. The capital Latin indices run 
••• = 0,1,2,3,5,6. 

The 4-dimensional Einstein equations are R^ffi} — ^iifi^R'''^^ = 0, i.e., the ordinary vacuum equations without no 
cosmological constant. The Greek indices a, /?,... = 0,1,2,3 refer to the ordinary four dimensions. The ansatz 
functions for the source are 

= -g,.uF{r), T,, = -g^jK{r), T,^ = . (169) 

The small Latin indices refer to the extra coordinates i, j = 5, 6. Finally the metric takes the form 

ds^ = (p^{r)T]^,ydx^dx'' - A(r)(dr2 + r^dd^) . (170) 

The metric of the ordinary 4-space, r\ap^ has the signature of (+, — , — , — ). The 4D and 2D "warp" factors, namely 
the ansatz functions 0(r) and A(r), depen d only on the extra radial coordinate, r. The Einstein field equations (jl68p 
in terms of the ansatz functions are |110l | 

3^ + 34P + ;^-;^ + ;^ = ^fF(r)-Ai , (itd 

r(f> 

(172) 
(173) 

where the prime ~ d/dr. These equations are the (a, a), (r, r), and {9, 0) components, respectively. Two of the three 
equations (|17ip - (|173p are independent, and they can be reduced to a set of two equations for (j){r), A(r). An analytic 
solution to these equations is found with the ansatz functions of the form 

^(^) - ^ A(r) . ^-^^ ^ -J , (174) 

where c and a are constants. The source functions are 
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where the constants /i = —^a and /2 = ^(a + 1) are determined by the 6D cosmological constant A, and the 
constant a, from the 4D warp function (j){r). 

A drawback of this solution is that the matter sources are put in by hand via the ansatz functions F{r) and K{r) 
rather than given on some realistic foundation of the field theory. Also from T^^, from (|169p and F{r) from (|175|) (as 
well as using the asymptotic values of F{r) at r = and r = oo) one sees that the energy density is negative on the 
brane and decreases to a negative value at r = oo. 

In the paper these solutions were generalized to the case of n > 2 extra dimensions. As with the 6D solution, 
the obtained solutions provide the universal, gravitational trapping mechanism of fields with various spins of from 
to 2. Considering the magnitude of the ansatz functions for the stress-energy which are necessary to realize this 
trapping solution one finds physically reasonable properties for a range of parameters n, A, e (brane width) and a. 
In addition, for n = 2 and rt = 3, it is possible to have the stress-energy tensor which approaches zero in the limit of 
r oo. For n = 4, either F{r) or K{r) goes to zero, but not both do. 



V. COSMOLOGICAL THICK BRANES 



In the previous sections, thick brane solutions which contain Minkowski or de Sitter four-dimensional geometry 
were considered. For such models, the bulk metric functions depend only on the radial bulk coordinate and solutions 
can be obtained, just by solving the ordinary differential equations. But if a thick brane which has a less symmetric, 
Fricdmann-Robertson- Walker (FRW) cosmological geometry, then in most of cases the metric functions also depend 
on time, and to find a solution, one has to solve par tial differential ec^uations. The cosmological thick brane solutions 
were not explored much, except for works e.g., [ll2|. 

In this section, several approaches to formulate the cosmological equations, in the presence of arbitrary matter 
inside the brane, are reviewed. Because of the subtlety about how to define and handle the effective four-dimensional 
quantities, there is no unique formulation. Thus, we review several representative approaches both in five- and 
multidimensional spacetimes. 



A. Codimension-one model 



1. Thin brane cosmology 



The interest in the localized matter distributions in the context of gravity has a long history. Because of the 
difficulty to treat a thick wall, in the early times it had been considered to idealize a wall as an infinitesimally 
thin object [liMil- 

Then, the interest had come again with the studies on cosmological phase transition s and 
formation of topological defects. Again, these defects were mainly assumed to be infinitesima lly t hin |ll6l . lll7l |. The 
formulation of the equations of motion of a singular domain wall was summarized by Israel |118| . In the context of 
the Lanczos-Darmois-Israel formalism, a thin shell is regarded as an idealized object with zero thickness. 

Recent developments of string and M-theory have motivated us to study cosmology in brane world. The concept of 
a domain wall was applied very frequently to treat a self-gravitating brane mathematically in the multidimensional 
general relativity. Cosmological equations on a thin brane distribution embedded into a five-dimensional bulk space- 
time is easily derived by inte grating th e five-dimensional Einstein equations with a (5-f unct i onal contribution into 
the energy-momentum tensor 119l 12Cll | or by employing the Israel junction conditions |l2ll . Il2 j |. The thin brane 



cosmology for a given bulk geometry is uniquely determined. As the simplest example, in an AdS-Schwarzschild bulk 
spacetime with the Z2 symmetry with respect to the brane, the effective Friedmann equation on the brane is given by 



a\ K 1 , 87rG4 / P\ M 
2 +^=3^^ + ^41+2^ ' 



+ — + ^ +-^, (176) 



where a and p are the scale factor and energy density of the brane, respectively, a represents the brane tension. The 
AdS curvature radius is given by i"^ = — 6/A, where A is a (negative) bulk cosmological constant. The constant M is 
related to the mass of a black hole sitting in the bulk. The four-dimensional Newton and cosmological constants are 
given by 

8,G. = f, A,.M„»-l = M.» + iA,. (177) 

In the effective Friedmann equation (jl76p . two corrections to the conventional one appear. The first correction is 
the term which is proportional to the square of energy density on the brane, which may be important in the early 
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universe. The second contribution is the so-caUed dark radiation which is proportional to the mass of the black hole 
M and a purely geometrical effect. 



2. Thickness of domain walls m general relativity 



In contrast to the case of a thin brane, the thickness brings new subtlety. Here, we see the brief history on the 
treatment of the thickness of a domain wall. Early attempts to formulate the thickness were mainly motivated in 
modeling the late time cosmological phase transition |l23l |. Silveira (l24| studied the dynamics of a spherical thick 
domain wall by defining an average of radius R, weighted by the energy density. It was also unclear whether the 
Israel conditions are applicable to the domain wall problem and the zero thickness limit of a thick domain wall is 
really equ i vale nt to a thin do main wall dominated by the Israel conditions. There were several works related to this 



125l - [l27 |. The work of jl28j | discussed a modification of the Israel thin shell equations in the Einstein-scalar 



issue 

theory and concluded that the effect of thickness at first approximation was effectively to reduce the energy density 
of the wall in comparison with the thin case, leading to a faster collapse of a spherical wall in vacuum. The work 
of [l29j | applied the expansion method of the wall action in terms of the small thickness parameter and showed that 
the effective action for a thick domain consists of the usual Nambu-Goto term and a contribution proportional to the 
induced Ricci scalar curvature. 



B. Thick brane cosmology 



Here, several representative approaches to discuss cosmology on a thick brane are reviewed. 

a. Averaging approach In the approach employed in Ref. |130{ . the effective brane quantities are obtained by 
integrating the five-dimensional ones over the brane thickness. The brane has a finite thickness between y = — yo/2 
and y = yo/2, where y is the proper coordinate along the extra dimension. For simplicity, it is assumed that the 
brane thickness is time-independent. 

There is ambiguity of the definition of the effective four-dimensional quantities. Here, the averaging prescription 
that one integrates a five-dimensional quantity Q{t, y) over the brane thickness 



'•yn/2 
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{Qit)) = - Q{t,y)dy, 

yO J-yo/2 



(178) 



is adopted. For convenience, the dimensionless quantities 
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(179) 



are introduced. The parameter r] represents the homogeneity of the matter distribution over the brane. By imposing 
Pr(±2/o/2) = 0, where Pr is the pressure along the extra direction, and by integrating the Einstein equations, the 
effective Friedmann equation is found to be 
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(180) 



where C is an integral constant. The constant C must satisfy the inequality that 

C(|<2o'(,-.)|)+»*('-|)- 

In the thin brane limit yo ^ 0, we find 



(181) 



(182) 
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which corresponds to the well-known thin brane cosmology with the quadratic dependence on the energy density of 
a brane. In the opposite limit ?/o ~^ c»i the generalized Friedmann equation reduces to 

(183) 

up to 0{y^^). In the case of a Minkowski bulk j/o ^ 0, the inequality yields C ^ and 




H^^'-^^pip), (184) 



which corresponds to the standard four-dimcnsional cosmology with effective Newton's constant k| :~ hi^rj/yo. In 
the homogeneous bulk 77 = 1, the result is in agreement with the Kaluza-Klein picture. Thus, the above equations 
smoothly interpolate between the Kaluza -Klei n and thin brane pictures of cosmology. The extension of this method 



to the DGP model was discussed in Ref. 131 



b. Two thin-shell approximation Another different approach based on the gluing of the regul ar thickened brane 
with two different spacetime geometries along two regular hypersurfaces was considered in [l32l Il33j . The purpose was 
to understand the dynamics of a localized matter distribution of any kind, confined within two different spacetimes or 
matter phases. It would be enable one to have any topology and any spacetime on each side of a thick wall or brane. 
By construction, the matching is regular and there is no singular surface whatsoever in this formulation. Therefore, 
the junction conditions for the induced metric extrinsic curvature tensors on the thick wall boundaries with the two 
embedding spacetimes can be applied as 

[/lahls. = 0' [^"afcls. - 0, (185) 

where (i=l,2) represents each reg ular boundar y loc ated between the core brane region and each bulk region. By 



applying the formalism developed in jl34| , in Ref. |135l | cosmology on the brane in an AdS-Sch bulk was investigated. 



It turns out that the first modified Friedmann equation becomes 

p 36 3 a* 

where the effective four-dimcnsional cosmological and gravitational constants can be read off 

A4 A w^A^ 4w{~A) 

T=6+^' 8-G4 = ^— . (187) 

2w is the proper separation between two regular walls. In contrast to the thin brane case, the linear dependence 
on the brane energy density is obtained even without the tension. Newton's constant is proportional to the brane 
thickness. Other than the energy density equation, one also obtains the pressure equation which depends on the bulk 
pressure in a nontrivial way. These two cosmological equations together with the energy c onser vation law determine 
the system completely. The time- variation of the brane thickness was discussed in Ref. |l36l |. It was shown that 
in the absence of the pressure along the extra dimension in the brane energy-momentum tensor the thickness of the 
brane decreases with time, while at a late time the negative transverse pressure can lead to an increase in the brane 
thickness. According to Eq. (|187p . this suggests the time variation of the constants G4 and A4. 



c. Quasi-static approximation The authors of Ref. [137| considered a thick codimension-one brane model in 
which the brane energy-momentum tensor includes the pressure component along the extra dimension. It is assumed 
that there is no matter outside the brane. By integrating the four-dimensional time and spatial components of 
the five-dimensional Einstein equations along the fifth dimension, the first-order transverse derivatives of the metric 
functions at the brane boundaries are related to the components of the brane energy-momentum tensor integrated 
over the brane. Here, the assumption that the derivatives parallel to the brane are negligible inside the thick brane 
in comparison with the transverse ones, is made. These matching conditions are plugged into the remaining extra- 
dimensional component of the Einstein equations and the cosmological acceleration equations become 
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where p = a + pm, P = — c + wpm and pr = —CTr + WrPm are energy density, pressure and bulk pressure in the brane, 
which are defined by the integration of the bulk stress tensor over the thick brane, respectively, a and (7^ are constant 
parts of the four-dimensional matter, i.e., tension and the bulk pressure, and {w,Wr) specify the equations of state 
of the time-dependent parts. The presence of the bulk pressure violates the energy conservation on the brane. The 
effective equation of state of the brane matter is given by 

W + Wr , 

W^eff = -■ (189) 

1 — 

An accelerating universe can be obtained for ifoff < —1/3. It imp lies that for example even for cosmic dust = 0, 
the cosmic acceleration can be realized as long as Wj. < —1/2 |l37| . 

This method was originally developed in the codimcnsion-two case and also generalized to the cases of higher 
codimensions, which will be reviewed in the next subsection. 



C. Higher-codimensional cases 

1. UV behavior of higher- codimensional branes 

In general, a 3-brane with more than two codimensions is a singular object because the gravitational potential 
sourced by the brane is divergent in approaching it. This is understood as follows. Imagine a codimension-n (n = 
1, 2, 3, • • • ) self-gravitating object. Gauss's law shows that the gravitational potential becomes V{r) cx l/r"~^, where 
r is the radial distance from the object. In the case of n = 1, which is a singular domain wall case, the jump across 
the object is finite and the matching conditions arc available. The case of n = 2, where the object is a string-like 
defect, is marginal and only the tensional matter can be put on the object and induces the deficit angle of the bulk. 
In the case of n > 3, for any kind of the localized matter, the potential becomes singular and a naked singularity is 
usually formed. In such a case, one needs some prescription, such as UV regularization, to put matter there. 

Several prescriptions for such divergence were suggested in the field theory approaches. In Refs. |l38l - ll4Cl| . an 
UV regularization of the self- forces of various fields (scalar, graviton and form fields), for extended objects with 
codimensions more than two, were discussed in the linearized approximation. The regularization was performed 
by replacing a singular delta functional source with a regular profile function (An IR cut-off is also required for a 
codimension two brane to cure the logarithmic divergence). For the gravitational force alone, the force is proportional 
to the extrinsic curvature vector, i.e., the trace of the second fundamental tensor. For the special case of a codimension- 
two brane, the force is actually zero and so the self-interactions cancel, which is a well-known result for cosmic strings 
in four dimensions, i.e., the tension determines the conical deficit angle and does not affect the local geometry both 
inside and outside the string. When a dilaton and a form field are included, certain combinations of couplings 
result in the vanishing total self-force. It w as a l so sh own that these regularized self-interactions can be expressed 
as renormalizations of the action. In Ref . |l4ll . Il42| , the classical renormalization schemes for a local divergence 
at a codimension-two (conical) brane were discussed, mainly for a scalar theory. It was shown that the logarithmic 
(classical) divergence of a bulk field in the thin brane limit can be renormalized into the brane localized mass and 
coupl ings, at the tree and one-loop levels, and the renormalization gro up e quations for them were obtained. In 



Ref. |14l| . the brane dynamics was not taken into account and in Ref. |142| it was done. The philosophy behind 
the procedure is that the low-energy physics on the brane should be independent of the high-energy regularization 
mechanism: 1) bulk fields evaluated away from the brane should not depend on the regularization mechanism and thus 
be finite in the thin-branc limit, 2) bulk fields evaluated on the brane itself would be sensitive to the regularization 
procedure and thus it would not be required that bulk fields evaluated on the brane are finite in the thin-branc limit, 
3) brane fields should have a well defined low-energy theory independent of the brane regularization. 



2. Cosmology on regularized higher-codimensional branes 

To discuss n onlinear g ravity and cosmology on a codim ension-two brane, various ways of regularization were derived, 
in e.g., Rcfs. [I43l4l50f . In the works Ref. |l43l Il44l |. the authors considered a way of regularizations to smooth 
out the matter over the brane with a finite thickness and investigated the effective cosmology in a class of six- 
dimensional models with stable and compact extra dimensions. The time-dependence of the brane matter was taken 
into consideration perturbatively around a static solution. By integrating over the brane, the effective cosmological 
equations are obtained and coincide with the conventional one at the low energy scales. 
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Another way of regularization of a codimension-two brane has been discussed in Refs. [Hi-liiii, such that a 
codimension-two brane was replaced with a regular bulk geometry, i.e., a capped region. Between the original bulk 
and the capped region, a ring-like 4-brane with the compact fifth dimension, wrapping around the axis of symmetry 
of the bulk, appears. This 4-brane is a codimension-one object and the Israel junction conditions can be applied. A 
motion of the brane realizes cosmology, but as long as the bulk is static, the effective cosmology is not the realistic 
one. An appropriate inclusion of the time-dependence of the bulk helps the recovery of it on the brane. 

The work |l50l | also discussed cosmology on a thick codimension-two brane. Apart from the previous approaches, 
the discussion of jl50l | was given in the general background and may be able to b e ap plied to various codimension-two 
br ane m odels, as long as several assumptions are satisfied. The authors of Ref. (l5l| extended the method employed 
in jl50l | to the case of an arbitrary number of codimensions. The discussion starts from the general metric ansatz that 



-N{t, rfdf + A{t, rfgijdx'dx^ + dr^ + a{t, rf-labdy^dy^ 



(190) 



where ^ab is the metric of (to — 1)— sphere. Following Ref. |l38l - ll40l |. regularization of the brane is performed by 
replacing the delta function with a regular profile function. As the profile function, here, it is assumed that all the 
brane matter is uniformly distributed over the brane region < r < e. 

Here several assumptions are going to be made: The important one is that the derivatives tangential to the brane 
are small enough in comparison with those normal to the brane. It is also assumed that A{t^ 0) sa A(t, e) = a{t) and 
the energy momentum tensor averaged over the brane takes the form of = diag( — p,p,p,p,Pr,Pb, ■ ■ ■ )i where p, 
p, Pr and ph are the energy density and pressures along the three-dimensional space, bulk radial direction and bulk 
angular directions, respectively. 

In the case of a codimension-two brane of m = 2. in Ref. the pressure equation is obtained as 
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Finally, in the codimension to > 3 cases, 
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where 
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and, Am = 2/ (fim-ia'" ^M™^^^ and = —2 — to. As in the five-dimensional case, the brane matter is decomposed 



into the constant (tension) part and timc-dcpcndcnt part as p — a ^ 
expanding the equation of motion up to the linear order pm, 



Pm and p = —(T + wp,n- In the low energy regime, 



(a^^) " (wi + 2CTV3 + 4w4(T) + (u;4 + cra;3)(l-3w)p™ + 0(p^„), 



(194) 



which has the form in the conventional cosmology. 

In the case of higher-codimensions that to > 2, under several assumptions such that the time evolution of the brane 
thickness is negligible, Pr is constant, and the energy flow into the bulk is negligible r^r-jg = 0, the standard energy 
conservation law is satisfied. The form is similar to the Randall-Sundrum one. But the effective gravitational constant 
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is proportional to which cxphcitly depends on the bulk pressures pi, and pr. However, in the higher codimensional 
case, there is the subtlety on the presence of a localized zero mode. 

In the case of a codimension-two brane m — 2, the difference from the cases of m > 2 comes from the term 
proportional to ^/—g\^_^/y/-^\^_f^- But by using the matching condition, as long as \pr\ and are small enough in 
comparison with \p\ and \p\, -\/— ff|^_g can be approximated to be unity. Under such an approximation, the 
cosmological behavior has no essential difference from the cases of m > 2, but there is no problem on the localization 
of a graviton zero mode. 



VI. S-BRANES 



Usually, a brane is considered as a time-like manifold. However, it is also possible to consider time-dependent 
solutions in multidimensional gravitationa l th eories and to interpret them as spacelike branes (S-branes). Such 
an approach was suggested in the paper |l52j . S-branes which are thick and regular everywhere are interesting 
cosmological solutions for solving the cosmological singularity problem. 

The authors of the paper |152| support their interest for investigations in this direction by the following arguments: 
SD-branes (S-branes with Dirichlet boundary conditions) may holographically reconstruct a time dimension just as 
in th e Ad S/CFT correspondence the D-brane field theory holographically reconstructs a spatial dimension. 

In [152| |. the physical situations giving the solutions of an S2-brane domain wall, an Sl-brane vortex, and a charged 
S-brane are considered. The arguments in favor of existence of S-branes in string theory are presented. 

Several solutions of the supergravity equations corresponding to S-branes with odd codimensions are obtained. The 
following SO-brane solution in the D=4 Einstein-Maxwell gravity is obtained. The starting action is 



{R - F^) . 



(195) 



The corresponding Einstein-Maxwell equations have the following solution 
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(196) 



where Q and tq are constants and is the unit metric on a 2-dimensional space with the constant negative 

curvature, respectively. The solution describes a charged SO-brane. 

An S5-brane in the D=ll supergravity is obtained as well. The starting action is 
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The solution of corresponding filed equations is 
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where q, to are constants and dH^ is the unit metric on a 4-dimensional space with the constant negative curvature, 
respectivel y. It was shown that there exists a singularity of this metric near the brane for a large t. 

In Ref. |153| S-branes containing a graviton, form field strength, Fj^j, and a dilaton scalar, 0, coupled to the 
form field with the coupling constant a, are discussed. In the Einstein frame, the action is given by 



S = J d^x^g {r - \d^cj^d^cp - ^f'^^Fy 
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The corresponding equations of motion are 
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A solution describing an S-brane is sought in the form 

ds^ = -e^^df + e^^idxl + ■■■+ dxl) + e^^ d^l, + e^^idyl + ■■■ + dyl_^). 
The hyperspace Sfc,o- is a maximally symmetric space with a constant scalar curvature 
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The corresponding field equations can be written in the following form 
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which describes S-branes. 

The S-brane solution obtained in Ref. |l54l | is based on the use of two strongly interacting scalar fields giving regular 
solutions over all spacetime. The solution is based on the fact that the potential of the scalar fields has two global 
and two local minima. As \t\ — > oo the scalar fields go to one of the local minima. This is the necessary condition 
for existence of such type of a solution. The numerical analysis showed that, apparently, there is no other solution 
which tend s toward one of the global minima. From the mathematical point of view, the obtained S-brane solution 
in Ref. |l54{ is very similar to the usual (time-like) brane solution considered in the section UlI Al 

The Lagrangian for gravitating phantom scalar fields 0, x is: 
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where the potential l^((/>, x) describing strong interaction between scalar fields is given by 
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where Vq is defined by initial conditions. The metric describing an S-brane is 
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After finding the corresponding equations for the gravitating scalar fields (j) and X: can see, from the physical 
point of view, that the metric describes the bounce of the universe at the moment t = 0. Corresponding equations 
(in the Planck units, i.e. at c = = G = 1) are 
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The numerical solution with the boundary conditions 

a(0) = ao, a(0) = 0, 

0(0) = 00, 0(0) = 0, 
X(0) = Xo, X(0) = 

gives the result presented in Figs. [TT] and [T^ It is necessary to note that, as well as for time-like branes described 
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FIG. 11: The graphs for the scalar fields (j)* (t) , x* (t) . Both 
functions tend asymptotically to the local minimum at <j)* (t) — 
ml,x''{t) = 0. 



FIG. 12: The Hubble parameter H — a/a tends asymptotically 
to the de Sitter stages: if ff < - the stage of contraction, if 
H > - the stage of expansion. 



in the section UlI Al regular solutions (presented in Figs. [TT] and [T^ exist only for definite values of the parameters 
nil = ml, 1712 = m2- In this sense, the solution of the equations (j216l) - (|219|) is a solution of a problem for two 
eigenvalues mi and m2 and two eigenfunctions 4'*{t) and X*{t)- 



VII. CONCLUSION 



In this paper, we gave a review of works devoted to study the thick brane solutions. Our main objective is to 
introduce the recent progress. As mentioned in the introduction (section I), there are the important motivations 
to consider the thick brane models: The first motivation comes from the interesting feature of a multidimensional 
topological defect that can trap the various fields inside it. The second one comes from the fact that along the 
developments of the (thin) brane cosmology, in particular in considering models of a spacetime with more than six 
dimensions, it turned out that even at the classical level a thin brane exhibits a singular behavior, due to the stronger 
self-gravity of the brane. To cure it, it would be crucial to take the effects of brane thickness into consideration. 
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We showed the previously known solutions and reviewed their basic properties. But as mentioned in the introdue- 
tion, in this article, little was referred to about the applications of these solutions to the high energy physics and the 
localizability of the fields with various spins, because they concern both the thin and thick brane models and should 
be reviewed in separate publications. We tried to introduce classifications of the thick branes on the basis of their 
division into some classes by combining these models on topological and physical properties . They are summarized 
below (see also tables I and II). 

Topological feature: The first important classification concerns whether a solution is topologically trivial or non- 
trivial. This is the division of the models in terms of topological properties. The detailed explanation about the 
concept of the topological triviality was given in the subsection II A. 

The presence of topologically nontrivial vacuum states is typical for the topologically nontrivial solutions (section 
II). Here there are two cases: 

1. codimension is equal to one 

2. codimensions are more than two. 

In the first case, there are kink-like solutions. Such solutions were found in the theory with scalar field(s) with a 
non-linear potential. In the second case, there are hedgehog-like solutions. Such a topological structure provides the 
stability of the solutions (at least at the classical level) that, of course, is a great advantage of such models. 

In the case of the topologically trivial solutions, there is only one vacuum state (see section III): a solution starts 
from a vacuum and returns to the same vacuum. The solutions were constructed by employing two strongly coupled 
scalar fields. A question concerning the stability still remains. The stability analysis should be done in case by case. 

Brane geometry: The geometry of the four-dimensional section on the brane is also an important point. The most 
works were devoted to search for a solution with a maximally symmetric four-dimensional section, i.e., Minkowski, de 
Sitter or anti-de Sitter(AdS) spacetime. This is because the system can be described by a set of the coupled ordinary 
differential equations with respect to the (radial) bulk coordinate. In particular, in the five-dimensional problem, the 
solution generation techniques as the superpotential method, discussed in the subsection II B, were developed. In a 
Minkowski or a de Sitter thick brane background, the trapping gravitons are possible, but in the AdS one, it may be 
impossible. The (stable) thick de Sitter brane solutions could be used for modeling inflation or dark energy. 

In terms of obtaining a realistic cosmological model, a solution with the cosmological four-dimensional section is 
more desired (see section V). They would be used for describing the evolution of a brane configuration changing with 
time. In this case, generically, the problem depends both on extra space-like coordinates and time. It is obvious that 
both the search and interpretation of solutions would be much more complicated. Therefore, the research in this field 
was performed by introducing the brane thickness by hand. Although there is always ambiguity on how to define and 
handle the thickness, several procedures were proposed, such as the averaging approach, the thin-shell approximation 
and the quasi-static approximation. In such models, generically, new effects due to the thickness of the brane are 
induced onto the effective Friedmann-like equations, which for instance could accelerate the Universe even though 
there is only the ordinary matter inside the brane. But the nature of these effects depends on the regularization 
procedure. 

Bulk geometry: In known solutions in five-dimensional models, the asymptotic bulk geometry can be either flat, 
de Sitter or AdS. In the model with a canonical scalar field, to obtain an asymptotically AdS thick brane solution, 
the fine-tuning condition Eq. (P7)) must be satisfied. In a multidimensional spacetime, the asymptotic geometry and 
regularity of the spacetime crucially depend on the boundary conditions and the choice of parameters. It is not only 
the case that the bulk spacetime is regular, but in some cases the bulk spacetime contains a naked singularity at a 
finite distance from the brane. 

Special models : The model with an unusual source function put in by hand provides another special kind of 
solutions (see section IV). In this case, the aim is to simplify a model as much as possible to obtain analytic solutions. 
It allows analyzing properties of the model in the most evident form. However, this advantage of the model is smoothed 
over by the fact that the source is put in by hand. 

The next special kind of solutions is the so-called S-brane model, in which it is assumed that the warp factor is 
some function of time (see section VI). It allows one to describe a cosmological evolution of multidimensional models 
which are both regular in a whole multidimensional space and do not have the initial cosmological singularity. 

Within the framework of multidimensional theories including string and Kaluza-Klein theories, the possibility of 
extra time-like dimensions may not be excluded. This may provide the disappearance of the cosmological constant 
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in the four-dimensional spacetime jl55l Il56l |. However, using the extra time-l ike dimen sions leads to considerable 
diffic ulties in fundamental properties, like appearance of tachyons and ghosts jl57l Il58| . and violation of causality 
[l59| . A model of a thick brane with extra time-like coordinate can be obtained, for example, with the use of two 
interacting scalar fields from section IIII Al In the case of a 5-dimensional model from IIII A 21 one can choose the 
metric in the form 

ds^ = 4)'^{r)r)^fi{x'')dx°'dxf^ - S dr^ , (224) 

where 6 = +1 corresponds to the problem with one time-like dimension considered in IIII Al and S = —1 refers to the 
problem with the additional time-like coordinate. In the latter case, the same thick brane solutions like in section 
IIII Al will exist if one takes the reverse sign for the potential (|157p . 



TABLE I: Classification of the thick brane models with scalar fields 



Class of brane models 


Matter sources 


Basic feature 


Topologically non-trivial 
thick branes 


Scalar field(s) with non- 
linear potentials 


The presence of topologically nontrivial vacuum 
states provide stability of the solutions at the clas- 
sical level. A question about the quantum stabil- 
ity requires further studies. 


Topologically trivial thick 
branes 


Two strongly interacting 
usual or phantom scalar 
fields 


A question about the classical and quantum sta- 
bility requires further studies. 



TABLE IL Special thick brane models 





Matter sources 


Basic feature 


Branes with unusual source 
functions 


Sources put in by hand via 
the special ansatz of func- 
tions 


The advantage is that it is possible to make use 
of analytical studies. The main drawback is that 
the source is put in by hand. 


S-branes 


Scalar fields, including fields 
with spontaneous symmetry 
breaking 


The multidimensional time-dependent solutions 
can be both topologically trivial and nontrivial 
ones. 


Thick branes with addi- 
tional time coordinates 


Two strongly interacting 
scalar fields 


They are not excluded but possess the consider- 
able difliculties such as appearance of tachyons, 
ghosts and violation of the causality. 



We shall close this review article after suggesting the possible perspective directions of investigations: 

1. To search new thick brane solutions, with the maximally symmetric four-dimensional section, in more generalized 
models of the field theory and the modified gravity. 

2. To search self-consistent solutions describing objects localized on the brane (monopoles, black holes and so on) 
depending both on extra space-like coordinates and 4-dimensional coordinates. 

3. To search the cosmological thick brane solutions which depend both on time and extra space-like coordinates. 
Having explicit solutions of such cosmological thick branes, the ambiguity mentioned above would be somewhat 
resolved. 
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